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Conversion  Factors, 

Non-Si  to  SI  Units  of  Measurements 


Non  SI  units  of  measurement  used  in  this  report  can  be  converted  to  SI  units  using 
the  conversion  factors  listed  below. 


Multiply 

By 

To  Obtain 

degrees  (angle) 

0.01745329 

radians 

feet 

0.3048 

metres  j 

indies 

0.0254 

metres  | 

inches  per  secord 
squared 

0.0254 

metres  per  second 
squared 

ksi  (kips  per  square  in.) 

6.894757 

megapascals 

pounds  (force) 

4.4484 

newtons 

1  pounds  (mass) 

0.4535924 

kilograms 

1  pounds  (force)  per 

1  square  inch 

0.006894757 

megapascals 

1  g  (standard  acceleration 

I  of  free  fall) 

9.80665 

metres  per  second 
squared 

1  Introduction 


Background 

Olmsted  Locks  and  Dam  (L&D)  is  one  of  the  largest  civil  works  projects 
undertaken  by  the  U.  S.  Army  Corps  of  Engineers  to  lead  the  modernization  of 
Corps  navigational  facilities  for  the  twenty-first  century.  Maintaining  a  robust 
navigational  infrastructure  to  facilitate  the  water  transportation  and  to  meet  the 
demand  for  the  ever-increasing  barge-traffic  flows  through  our  nation’s  inland 
waterways  is  vital  to  our  economy.  Transportation  of  bulk  commodities  through 
our  vast  inland  waterways  provides  the  most  economic  mode  of  conveyance  and 
it  conserves  natural  energy  resources.  In  this  complex  global  economy, 
increasing  productivity  by  making  efficient  and  effective  use  of  resources  has 
become  increasingly  important.  In  this  regard,  the  Corps  continuing  effort  to 
improve  the  navigation  facilities  using  the  latest  technology  is  exemplary. 

An  experimental  and  analytical  investigation  was  initiated  by  the  U.  S. 
Army  Engineer  (USAE)  Division,  Ohio  River  (ORD)  and  the  USAE  District, 
Louisville  (LRL)  to  study  the  flow-induced  characteristics  of  the  wicket  gates 
proposed  for  the  Olmsted  navigation  dam.  Several  research  and  development 
phases  were  initiated  and  coordinated  by  the  LRL,  to  accomplish  this 
navigational  project  on  the  Ohio  River.  The  principal  focus  of  this  research 
scheme  was  to  determine  the  most  appropriate  design  of  wickets  for  the  new 
Olmsted  L&D  project.  Scale  models  based  on  hydraulic  and  structural 
similitude  relationships  have  been  studied  to  understand  the  flow-induced 
dynamic  characteristics  of  Olmsted  wicket  installed  in  the  Olmsted  Prototype 
Wicket  Dam,  Smithland  facility.  These  models  have  provided  researchers  with 
necessary  tools  for  the  structural  performance  evaluation  of  hydraulically 
actuated  wickets  (March  and  Elder  1992,  Seda-Sanabria  1994,  and  Chowdhury, 
Hall,  and  Pesantes  1997).  Such  evaluation  requires  the  development  of  complex 
analytical  models  capable  of  properly  simulating  the  interaction  between  the 
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structure  and  the  fluid  flow- fields.  For  this  purpose,  the  ABAQUS/Standard  FE 
module  was  used  to  develop  a  numerical  model  of  the  Olmsted  wicket. 

The  inability  of  the  physical  models  to  address  broad  boundaries  of 
enquiry  for  numerous  operating  conditions  demanded  the  development  of  an 
analytical  model  at  WES.  This  report  presents  the  step-by-step  updating 
procedures  used  for  developing  an  FE  model  based  on  dry  and  wet  experimental 
modal  analyses  and  flow-induced  dynamic  experimental  results.  An  updated 
model  would  enable  the  researchers  to  study  the  effect  of  future  modifications,  if 
any,  in  obtaining  an  efficient  and  optimum  design  for  the  proposed  wicket.  This 
model  can  be  used  to  measure  the  sensitivity  of  different  parameters  influencing 
the  dynamics  of  the  wicket  and  to  evaluate  the  structural  dynamic  performance 
during  critical  operation. 

Olmsted  Navigational  Pass 

The  Olmsted  L&D  project  is  to  be  constructed  on  the  Ohio  River,  16 
miles  upstream  from  the  confluence  of  the  Ohio  and  the  Mississippi  Rivers. 

This  new  L&D  project  will  replace  L&D  52  and  53  on  the  lower  Ohio  River. 

This  facility  will  facilitate  navigation  at  Paducah,  Kentucky,  with  a  capacity 
sufficient  to  meet  projected  demands  for  tow  traffic  through  the  year  2025. 
Construction  of  the  Olmsted  L&D  Project  was  authorized  by  the  United  States 
Congress  on  17  November  1988,  by  the  passage  of  the  Water  Resources 
Development  Act  of  1988  (Public  Law  100-676).  The  cost  of  this  project  is 
being  equally  shared  with  the  navigation  industry.  Tariffs  paid  by  the 
navigational  traffics  on  diesel  fuel  are  used  to  form  an  Inland  Waterways  Trust 
Fund  which  will  provide  50  percent  of  the  project  cost.  Estimated  total  project 
cost  is  over  $1  billion. 

Wicket  gates  have  been  used  as  part  of  lock-and-dam  systems  by  the 
Corps  of  Engineers  since  1909  (Soast  1994).  Originally,  220  remotely-operated 
hydraulic  actuators  were  proposed  to  raise  and  lower  the  wickets.  However,  the 
current  plan  calls  for  a  1 ,400-ft  section  with  boat-operated  steel  gates  in 
combination  with  six  tainter  gates  and  a  fixed  weir.  At  10-ft  wide,  26  -ft  long. 
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and  with  a  design  lift  of  21 -ft,  these  wickets  will  be  the  largest  hydraulically- 
actuated  wickets  in  the  world.  The  LRL  gate  design  team  has  designed  the 
prototype  model  for  the  Olmsted  wickets  built  in  the  Smithland  Dam  on  the 
Ohio  River.  The  project  will  include  twin  1,200-ft  by  1 10-ft  lock  chambers  with 
a  design  lift  of  21 -ft  and  a  1,400-ft  long  navigation  pass  and  660  ft  of  tainter 
gates.  In  the  raised  position,  the  wicket  gates  will  be  used  to  regulate  the  level 
of  the  navigable  pool.  In  the  lowered  position,  the  gates  will  provide  a  navigable 
pass  (Lance  1992).  The  remaining  140-foot  dam  will  include  a  fixed- weir 
section  which  will  tie  into  the  Kentucky  shore. 

The  Olmsted  navigation  pass  is  one  of  the  busiest  barge  traffic  route  in 
the  nation.  The  new  L&D  project  will  be  constructed  near  the  community  of 
Olmsted,  IL,  at  Ohio  River  Mile  964.4.  Figure  1  shows  the  relative  location  of 
the  Olmsted  Dam,  the  navigational  routes,  and  the  old  L&D  52  and  53.  As 
shown  in  the  figure,  this  reach  of  the  Ohio  River  is  particularly  strategic  in  that  it 
provides  a  connection  between  the  Ohio,  Tennessee,  Cumberland,  and 
Mississippi  Rivers.  The  area  has  been  described  as  the  "hub"  of  the  Ohio  and 
Mississippi  Rivers  waterways  system.  Barge  traffic  moving  between  the 
Mississippi  River  system  and  the  Ohio,  Tennessee,  and  Cumberland  Rivers  must 
pass  through  this  stretch  of  river.  More  tonnage  passes  this  point  than  any  other 
place  in  the  inland  navigation  system.  In  1991  alone,  100  million  tons  of  goods 
were  shipped  through  this  junction  of  the  Ohio  River.  Clearly,  this  is  a  critical 
reach  of  water  from  a  commercial  navigation  perspective. 


Figure  1 .  Olmsted  L&D  location 
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The  original  dams  consisted  of  wooden  wickets  that  were  manually 
raised  to  hold  back  water  during  periods  of  low  flow  and  dropped  to  the 
river  bottom  during  high  water.  A  view  across  the  Ohio  River  at  L&D  No. 
52,  showing  the  operational  condition  of  the  wooden  wickets  is  displayed 
in  Figure  2.  Built  in  1929,  these  aged  L&D  were  renovated  to  meet  traffic 
demands  in  1970's  by  adding  temporary,  1200-ft  new  lock  chamber.  The 
present  facility,  however,  is  inadequate  to  handle  the  current  traffic 
volume  without  significant  delays.  The  continuing  growth  in  barge  traffic 
on  the  Ohio  River  and,  more  importantly,  deterioration  of  these  aging 
structures  demand  a  replacement  of  the  L&D  in  this  strategic  hub  of  the 
waterways. 


Figure  2.  Locks  and  dam  #52 

The  current  operating  proposal  calls  for  the  wicket  gates  to  be 
raised  in  groups  of  five  with  five-gate  gaps  until  all  110  gates  are  raised. 
Flow  over  the  top  of  any  gate  is  not  anticipated  except  during  its  raising  or 
lowering  operation.  The  five-gate  gaps  will  then  be  closed  by  raising  two 
gates  simultaneously  or  one  gate  at  a  time.  In  each  case,  the  end  gate  or 
gates  will  be  moved.  To  maintain  a  uniform  flow  pattern  in  the  river,  as 
many  as  44  gates  will  be  simultaneously  and  continuously  operated  in  a 
manner  that  subjects  the  gates  to  uneven  hydraulic  loading  that  may 
produce  flow-induced  gate  vibrations.  In  addition,  the  raising  or  lowering 
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operations  could  also  induce  gate  vibrations.  A  large  amount  of  gate 
manipulation  is  expected  during  low-flow  periods  due  to  powerhouse- 
controlled  inflows  into  Olmsted  Reservoir  from  the  Tennessee  and 
Cumberland  Rivers  (Elder  1992). 

An  updated  FE  model  capable  of  simulating  the  flow-induced 
vibrational  response  will  enable  the  LRL  to  evaluate  the  gate  performance 
for  a  variety  of  channel  and  flow  configurations.  The  analytical  capability 
saves  time  and  cost  by  reducing  experimental  simulation  of  numerous  test 
cases.  Most  importantly  the  long  term  performance  of  operational  wickets 
can  be  evaluated  with  this  FE  model  -without  much  additional  cost.  Gate 
component  performances  and  the  impact  of  individual  components  on  the 
overall  survivability  of  the  wicket  would  enable  the  district  to  measure  the 
reliability  of  the  system  as  well.  The  findings  obtained  from  these  analyses 
will  help  in  the  support  of  an  effective  design  procedure  for  hydraulically 
lifted  wickets. 

Physical  Model  Studies 

The  Olmsted  Wicket  Model  studies  began  in  January  1 990  when 
the  LRL  tasked  the  WES  with  the  development  of  a  1  ;25-scale  curved  gate 
model  of  the  Olmsted  L&D.  Since  then,  two  additional  models  were 
developed;  1 :25-  and  1 :5-scale  flat  gates.  These  physical  models  have 
been  used  to  investigate  the  flow-induced  vibrational  response  of  Olmsted 
Wicket,  and  have  also  helped  researchers  to  better  understand  the 
complexities  involved  in  the  fluid-structure  interaction  problem  affecting 
them.  The  Olmsted  Wicket  Models  studies  also  included  a  prototype  test 
facility  in  the  Smithland  Dam  on  the  Ohio  River  (Chowdhury,  Ross,  and 
Hall  1 997).  This  facility  examined  the  physical  performances  of  various 
mechanical/electrical/hydraulic  systems  proposed  for  the  Olmsted  Dam. 
Although  the  prototype  facility  and  the  1 :25-scale  models  provided  the 
necessary  information  about  wicket  gate  vibrational  behavior,  both  of 
these  models  had  limitations.  The  prototype  model  was  limited  in  the 
flow  conditions  that  it  can  simulate,  and  the  1 :25-scale  models  were  not 
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true  similitude  structural  models.  The  1 :25-scale  models  and  the  prototype 
test  facility  and  their  limitations  are  further  discussed  by  Chowdhury,  Hall 
and  Davis  (technical  report  (TR)  in  preparation). 

Considering  the  limitations  involved  in  the  construction  of  a 
prototype  facility  (i.e.,  operational  inconveniences,  excessive  modification 
costs,  etc.),  and  the  inadequacy  of  the  1 :25-scale  flat  gate  model  in 
representing  the  structural  similarity  of  the  prototype  wicket,  it  was 
imperative  to  study  the  behavior  of  an  intermediate-scale  wicket  in  an 
easily  controlled  environment.  Thus,  a  moderate-scale  of  1 :5  of  the 
prototype  geometry  was  considered  for  fiirther  investigation.  The  1 :5- 
scale  flat  gate  model  was  capable  of  simulating  the  structural  as  well  as  the 
hydraulic  responses  of  the  prototype  gate.  The  1 : 5 -scale  flat  gate  model, 
shown  in  Figure  3,  was  geometrically  and  dynamically  a  similitude  model. 
In  contrast  to  the  previous  1 :25-scale  model,  its  geometrical  characteristics 
simulated  those  found  on  the  prototype  structure. 


Figure  3.  Olmsted  1 :5-scale  model 

Experimental  results  obtained  from  the  1 : 5 -scale  model  have  been 
used  for  updating  the  FE  model  presented  in  this  report.  These  included 
the  static  and  dynamic  parameters  indicative  of  the  actual  performance  of 
the  operating  wicket. 
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Objectives 


The  objectives  of  the  l:5-scale  flat  gate  model  investigation  in  support  of 
the  Olmsted  Model  studies  program  were  to: 

a.  Develop  an  analytical  finite-element  (FE)  model  for  predicting  the 
flow-induced  vibrational  response  of  Olmsted  wicket  gates  using 
analytical-  experimental  correlation. 

b.  Validate  design  of  the  wickets  by  evaluating  the  FE  model. 

c.  Investigate  the  gate  behavior  using  the  FE  model. 

Scope 

An  experimental  and  analytical  investigation  was  initiated  to  study  the 
flow-induced  characteristics  of  the  wickets  proposed  for  the  Olmsted  navigation 
dam.  Scale  models  based  on  hydraulic  and  structural  similitude  relationships 
have  been  studied  to  understand  the  flow-induced  dynamic  characteristics  of 
Olmsted  wicket  installed  in  the  Smithland  prototype  facility.  These  models  have 
provided  sufficient  experimental  data  for  developing  analytical  models  capable 
of  properly  simulating  the  interaction  between  the  structure  and  the  fluid  flow- 
fields.  For  this  purpose,  the  ABAQU S/Standard  FE  module  was  used  to  develop 
a  numerical  model  of  the  Olmsted  wicket.  The  FE  model  was  validated  using 
results  obtained  from  both  dry  and  wet  experimental  modal  analyses,  and  static 
and  flow-induced  dynamic  experiments  (Chowdhury,  Hall,  and  Davis  TR  in 
preparation).  Geometric  dimensions,  physical  properties,  and  experimental 
modal  data  were  correlated  with  the  analytical  model  to  ensure  the  correct 
distribution  of  mass  and  stiffness  of  the  model.  Wet  modal  parameters  were  used 
to  define  the  added  mass  distribution  of  the  operational  gate.  Flow-induced 
operational  test  pressures  on  the  gate  surfaces  were  used  to  define  the  loading 
function  and  its  history  for  the  FE  model.  The  final  product  of  these  studies  is 
an  updated  FE  model  based  on  these  correlations  that  adequately  predicts  the 
structural  response  of  the  Olmsted  wickets  for  various  operational  and  flow 
configurations.  This  report  is  a  review  of  the  step-by  -step  updating  process 
used  to  develop  the  final  FE  model.  It  also  describes  the  uses  of  the  model  for 
evaluating  the  dynamic  performance  of  the  operational  wicket. 
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2  Finite-Element  (FE)  Model 
Development  Phases 


Background 

The  ABAQUS  general  purpose  FE  code  had  been  used  to  develop  a 
numerical  model  for  the  Olmsted  wicket  structure  presented  in  this  report.  The 
general  purpose  of  the  model  was  to  simulate  the  flow-induced  structural  motion 
of  the  wicket  due  to  various  operating  conditions.  This  required  a 
comprehensive  modeling  of  the  fluid-structure  interaction  of  the  system  using 
experimental  results  so  that  a  Lagrengian  based  FE  code  could  be  used  for 
dynamic  response  prediction  of  an  operational  gate.  Experimentally  determined 
Lagrangian  displacement  based  fluid  elements  could  be  readily  incorporated  into 
a  general  purpose  computer  code  for  structural  analysis  (Calayir  and  Dumanoglu 
1993).  Experimental  results  have  been  extensively  used  by  the  researchers  to 
compute  a  virtual  added  mass  representing  the  fluid-structure  interaction  of 
major  navigational  structures  (Daniell  and  Taylor  1994;  Ishii  and  Knisely  1992; 
and  Thang  1982).  An  equivalent  added  mass  is  computed  such  that  the  free 
vibrational  natural  characteristics  of  the  wet  system  are  preserved  in  the  model. 

In  the  Lagrangian  approach,  the  displacement  variables  in  both  the  fluid 
and  the  structure  are  identical.  Thus  by  simulating  the  water  as  a  mass  added  to 
the  dry  structure  the  compatibility  and  equilibrium  are  automatically  satisfied  at 
the  nodes  along  the  interfaces  between  the  fluid  and  the  structures. 

Experimental  data  for  the  Olmsted  1 :5-scale  model  were  used  to  obtain  the 
added  mass  of  the  wicket  for  different  operating  conditions.  Experimentally 
recorded  known  pressure  fields  provided  the  input  forcing  functions  for  the 
operational  wicket.  This  completed  the  modeling  of  the  operational  wicket  in 
three  major  steps.  In  the  first  step,  the  geometric  model  was  developed  to  match 
the  mass  and  stiffness  of  the  dry  wicket.  In  the  second  step,  added  masses  were 
distributed  in  the  geometry  to  simulate  the  wet  natural  characteristics  of  the 
operational  gate,  and  finally,  the  forcing  functions  were  applied  to  the  FE  model 
to  measure  the  dynamic  response  of  the  operating  gate. 

Static  and  dynamic  modal  updating  were  used  to  calibrate  the  mass  and 
stiffness  of  the  FE  model  with  those  of  the  experimental  wicket.  Experimental 
modal  data  have  been  extensively  used  by  researchers  to  correct  the  FE  model 
for  vibration  analysis  (Okuraa  1993;  Niedbal  1985;  Touhei  and  Ohmachi  1994; 

Dascotte  1991;  Petrick  1993;  Brillhart  et  al.  1991;  and  Ewins  1984).  All  of  these 
researchers  except  Touhei  and  Ohmachi  (1994)  used  experimental  data  for 
updating  their  FE  model.  Touhi  and  Ohmachi  (1994)  used  a  finite  element- 
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boundary  element  method  in  the  time-domain  to  determine  the  modal 
characteristics  of  dam-foundation  system.  Updating  techniques  used  by  these 
researchers  included  the  modal  assurance  criteria  (MAC)  and  cross¬ 
orthogonality  check  (COC).  MAC  is  used  to  check  the  consistency  of  the 
extracted  modal  parameters  by  measuring  the  correlation  between  analytical  and 
experimental  modal  vectors.  The  COC  determines  the  orthogonality  of 
analytical  and  experimental  modal  parameters  based  on  the  property  of 
orthogonality  with  respect  to  the  analytical  mass  matrix. 

An  updated  FE  model  is  most  often  an  approximate  analysis  tool  which 
tends  to  minimize  the  deviation  between  the  analytical  and  experimental 
responses.  In  reality,  this  deviation  widens  due  to  the  uncertainty  and 
randomness  of  the  geometry  and  material  properties  of  the  structure,  and  loading 
and  boundary  conditions  (Chang  1993).  This  is  due  to  the  fact  that  in  most  FE 
analysis,  the  geometric  and  material  properties,  and  loading  conditions  are 
treated  as  deterministic  despite  the  fact  that  they  are  variant.  Difficulty  of 
modeling  of  boundary  conditions  becomes  one  of  the  most  troubled  area  in  most 
FE  analysis.  Such  an  inherent  limitation  prompted  researchers  to  study  the 
sensitivity  of  major  design  parameters  to  determine  the  structural  modification 
effects  in  the  analysis  for  many  decades  (Lin  Du  and  Ong  1993;  and  Bretl  1992). 
A  sensitivity  based  study  could  correctly  identify  the  controlling  parameters 
which  could  be  used  with  greater  confidence  to  measure  the  bounds  of 
vibrational  responses.  In  this  research,  a  sensitivity  study  was  conducted  to 
determine  the  major  influential  parameters  and  the  bounds  of  wicket  responses 
were  computed  using  the  most  probable  limits  of  boundary  and  loading 
conditions. 


Physical  Model 

Model  description 

Test  gate  #6,  shown  in  Figure  3,  was  used  to  extract  the  experimental 
data  for  the  FE  model  updating.  General  dimensions  and  the  structural 
components  used  for  fabricating  the  wicket  are  presented  in  Appendix  A.  A 
hydraulic  lifting  cylinder  was  used  to  raise  and  lower  the  gate,  and  the  prop  rod 
supported  the  gate  when  the  dam  was  in  service.  Lugs  mounted  on  the  sill 
bearing  frame  provided  base-hinge  mounting  connections  for  the  gate.  Base 
hinges  restrained  the  translational  movement  of  the  gate  with  free  rotational 
capabilities.  The  hydraulic  system  consisted  of  a  primary  and  an  alignment 
cylinder.  A  cup  on  the  end  of  the  primary  cylinder  piston  rod  provided  a  non- 
rigid  bearing  connection  (roller)  to  the  gate  through  a  ball  mounted  to  the  back 
of  the  gate.  This  roller  attachment  allowed  the  retraction  of  the  piston  rod  when 
the  gate  was  being  supported  by  the  prop  rod. 

One  end  of  the  prop  rod  was  pin  connected  to  the  gate  and  the  other  end 
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rested  on  the  hurter  as  shown  in  Figure  4.  The  bottom  end  of  the  prop  blade 
provided  "no-tension"  mechanism  such  that  the  prop  blade  was  restrained  only 
against  sliding  away  from  the  gate.  Top  end  of  the  clevis  connection  to  the  ball 
mounting  plate  provided  no  rotational  restraint.  As  a  result,  the  top  end  of  the 
prop  rod  randomly  shifted  its  position  during  the  operation  of  the  gate.  Such 
uncertain  behaviors  of  the  prop  ends  introduced  nonlinearity  into  the  system 
along  with  time-invariant  modal  responses  (Chowdhury,  Hall,  and  Davis  in 
preparation).  There  are  two  one-way  tracks  in  the  hurter  (see  appendix  A  and 
Fig  4b)  which  provide  the  sliding  paths  for  the  prop  rod  end-blade  during  raising 
and  lowering  cycles. 


a.  Ball-cup  and  prop-pin  connections  b.  In-place  prop-blade  on  the  hurter 

recess 


Figure  4.  Views  of  connection  details 

During  a  raising  cycle,  the  primary  cylinder  was  extended,  thereby 
rotating  the  gate  through  an  arc  from  -3  to  68  deg,  then  stopped  and  retracted  and 
the  load  from  the  gate  was  transferred  to  the  prop  rod  at  the  65-deg  position 
where  the  prop  rod  locks  in  the  hurter.  During  a  lowering  cycle,  the  primary 
cylinder  was  extended  while  the  alignment  cylinder  was  used  to  rotate  the 
primary  cylinder  to  match  the  cylinder  cup  with  the  ball  mounted  to  the  back  of 
the  gate,  the  gate  was  then  rotated  from  65  to  70  deg  where  the  prop  rod  was 
disengaged  and  the  primary  cylinder  was  retracted  until  the  gate  came  to  rest  at  a 
-3  deg  position. 


Physical  and  material  properties 

Prototype  physical  dimensions  shown  in  appendix  A  were  used  to 
construct  the  FE  model  presented  in  this  report.  Table  1  shows  the  elastic 
properties  of  major  structural  components  used  for  the  scaled  physical  model.  In 
this  table,  column  two  presents  the  design-drawing  reference  sheet  numbers 
corresponding  to  each  component  mentioned  in  the  first  column.  Tensile 
strength  experiments  were  conducted  for  each  batch  of  materials  procured  for 
the  scaled  physical  model.  Tensile  strength  results  for  each  component  are  also 
shown  in  the  table.  All  material  experiments  were  conducted  in  the  Structures 
Laboratory,  SL,  using  the  ASTM  E8  Test  Methods  (ASTM  1995).  As  seen  in 
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the  table,  material  properties  differed  from  components  to  components  and  even 
within  the  group  the  material  properties  varied  randomly.  An  average  of  three 
specimen  results  were  used  to  obtain  the  elastic  properties  presented  in  this 
table. 


Table  1 

Material  Properties  for  the  1:5-Scale  Wicket  Model 

Component 

Description 

Reference 
Sheet  # 

Material  Type 

Young’s 
Modulus,  ksi 

Yield  Strength,  psi 

Skin  Plate 

S-24,  S-25 

A  606 

40,974 

64,521 

Beam 

S-24,  S-25 

A  588  Grade  50 

34,287 

55,956 

Channel 

S-24 

A  588  Grade  50 

32,213 

55,956 

Stiffener 

Plates 

S-24 

A  514  Grade  70 

42,083 

61,216 

Web 

Stiffener 

S-24 

A  514  Grade  70 

32,480 

65,469 

Gate  Frame 
Hinge  pins 

M-8 

A  276  Stainless 
Steel 

30,850 

50,676 

Gate  Frame 
Hinge 

M-8 

A  588  Grade  50 

35,658 

52,128 

Ball 

Mounting 

Plate 

M-15 

A  276  Stainless 
Steel 

30,850 

50,676 

Prop  Rod 

M-18 

A  588  Grade  50 

30,737 

64,940 

Using  the  Archimedes  principle  of  buoyancy,  the  density  of  the 
wicket  was  measured  by  submerging  the  gate  in  a  water-filled  tank  and 
measuring  the  volume  of  water  displaced.  The  density  for  the  wicket  gate  shown 
in  Figure  5  was  421.98  Ib/ft^.  The  center  of  gravity  (CG)  of  the  wicket  was  also 
determined  by  balancing  the  gate  over  a  roller  as  shown  in  Figure  5.  The  CG  for 
the  wicket  was  24.875  in.  away  from  the  base  of  the  gate. 
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Figure  5.  Equilibrium  experiment  on  the  physical  wicket  model 


Preliminary  FE  Model 

The  preliminary  FE  model  of  the  test  wicket  was  developed  from 
construction  drawings  and  observation  of  the  1 :5  scale  structure.  This  initial 
model  represented  the  wicket  geometry  independent  of  its  supporting  devices  as 
shown  in  Figure  6.  Detailed  modeling  of  the  supporting  mechanisms  was 
accomplished  in  a  subsequent  study  presented  in  the  later  part  of  this  report. 

This  FE  model  included  1 1,763  nodes  with  a  total  of  3,179  elements.  Most  of 
the  structural  components  of  the  wicket  have  been  modeled  using  eight-noded 
shell  elements  (S8R5).  These  are  second-order  quadrilateral  shells  with  five 
degrees-of-freedom  (DOF)  per  node.  Six-noded  triangular  shell  elements 
(STRI65)  were  used  as  transitional  elements  to  provide  connectivity  of 
dissimilar  intersecting  elements  in  the  model;  these  elements  also  have  five  DOF 
per  node.  A  few  of  the  four-noded  quadrilateral  shell  elements  (S4R5)  were  also 
used  as  a  transitional  element  to  provide  continuity  of  the  nodal  pattern.  Two- 
noded  linear  beam  elements  (B31)  were  used  to  model  steel  rod  supporting 
hanger  and  the  instrumentation  cable  routing  pipes.  Lumped  mass  elements 
were  used  to  model  the  mass  of  both  bearing  brackets  and  the  connecting 
supports  of  the  hydraulic  cylinder  (see  Appendix  A  and  Fig  4). 

Gate  hinges  were  modeled  as  pinned  supports  and  the  prop  rod  was 
simulated  by  a  spring  applied  to  the  gate  bottom  surface.  Nodes  connecting  the 
gate  to  the  hinge  assemblies  were  restrained  at  both  hinges.  For  these  nodes, 
only  the  translational  DOF,  about  the  Z  direction  were  restrained.  Two  nodes  at 
the  hinges  were  also  used  to  restrain  the  translational  movement  of  the  hinge 
assemblies  about  the  X-and  Y-directions.  These  allowed  the  rotational  freedom 
of  the  hinge  assembly  and  restraint  about  the  Z-direction.  For  the  pin-connected 
hoisting  mechanism,  translational  DOF  about  Z-direction  were  restrained  for 
two  nodes  on  the  ball-mounting  plate  at  the  pillow-block  locations. 
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The  spring  coefficient  was  calculated  based  on  the  dimensions  of  the 
prop  rod  and  the  elastic  modulus  of  steel  and  was  1,822,974  Ibf/in.  Since  the 
support  elements,  which  included  the  gate  and  prop  rod  hinges  and  the  pillow 
block  and  ball  assembly,  contributed  significantly  to  the  mass  of  the  gate,  the 
mass  contribution  of  these  elements  was  calculated  and  applied  to  the  FE  model 
at  the  appropriate  nodes.  The  initial  FE  model  is  shown  in  Figure  6. 


Verification  of  Preliminary  Model 

Initial  verification  of  the  FE  model  included  physical  properties 
comparison  of  both  analytical  and  experimental  results.  The  total  mass  and  the 
center  of  gravity  for  both  analytical  and  experimental  models  are  compared  in 
Table  2.  A  close  agreement  between  the  predicted  and  actual  parameters 
indicates  that  the  mass  distribution  for  both  cases  were  identical. 

Further  correlation  of  the  preliminary  FE  model  consisted  of  comparing 
results  of  modal  dynamic  response  and  static  deflection  analysis  with 
experimental  results.  Experimentally  measured  modal  parameters  for  the  prop- 
supported  gate  were  compared  with  that  of  the  analytically  computed  data. 
Experimental  modal  testing  using  the  state-of-the-art  Scanning  Laser  Doppler 
Vibrometer  (SLDV)  system  was  conducted  to  extract  the  dry  and  wet  vibrational 
characteristics  of  the  1 :5-scale  physical  wicket  model.  Modal  testing  was 
conducted  by  exciting  the  wicket  gate  with  a  shaker  and  scanning  the 
downstream  surface  of  the  gate  using  the  SLDV  system.  Detailed  information 
on  data  acquisition  and  analysis  techniques  are  summarized  in  Chowdhury,  Hall, 
and  Davis  (in  preparation). 


Table  2 

Physical  Properties  for  the  Experiment  and  Analytical  Models 

Descriptions 

Experimental 

FE  Model 

Percent  Error 

Center  of  mass 
measured  from 
the  base  of  the 
gate,  y-dir 

24.84  in. 

24.57  in. 

1.09 

Total  mass 

182.75  lb 

182.38  1b 

0.20 
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Figure  6.  Finite  element  grid 


A  good  correlation  of  the  free-vibrational  characteristics  of  both 
analytical  and  experimental  models  would  indicate  a  dynamic  similarity  of  the 
models  which  in  turn  would  validate  that  the  mass  and  stiffness  of  the  actual 
gate  are  properly  modeled.  In  addition  to  the  dynamic  similarity  verification  of 
the  models,  a  static  load-deflection  test  was  also  used  to  calibrate  the  static 
stiffness  of  the  wicket. 

In  the  static  tests,  the  gate  was  placed  in  the  horizontal  position  resting 
on  the  supporting  hydraulic  cylinder.  Various  loadings  were  applied  to  the  top 
surface  of  the  gate  and  displacements  near  the  hinges,  at  the  center  and  near  the 
top  of  the  gate  were  measured  for  each  loading  condition.  In  the  dynamic  flow 
experiments,  a  single  gate  near  the  center  of  the  model  was  instrumented  with  a 
number  of  displacement  gauges,  pressure  gauges,  load  cells,  and  accelerometers, 
Pressures  on  the  gate  surfaces,  reaction  forces,  and  gate  accelerations  were 
measured  for  different  flow  and  channel  configurations.  Experimental  data 
recording  and  analysis  are  further  discussed  by  Chowdhury,  Hall,  and  Davis  (in 
preparation).  Flow-  induced  operating  responses  for  the  scaled  physical  model 
were  measured  for  different  gate  configurations  with  appropriate  head  of  pool 
elevations. 
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A  linear-elastic  technique  for  extraction  of  eigenvalue  data  for  the 
structural  system  was  used  to  compute  the  natural  frequencies  of  the  FE  model.  A 
detailed  formulation  for  the  ABAQUS  eigenvalue  and  other  dynamic  response 
analyses  methodologies  is  presented  in  the  literature  (HKS  1994).  The  response 
spectrum  analysis  of  the  initial  FE  gate  model  for  the  “dry”  condition  extracted  10 
modes  of  vibration.  Table  3  presents  the  analytically  computed  modal  frequencies 
for  the  wicket  model  for  two  different  types  of  boundary  conditions.  Different 
boundary  conditions  were  examined  to  measure  the  sensitivity  of  supporting 
mechanisms  on  the  natural  dynamics  of  the  wicket.  A  free-free  supporting 
condition  represented  the  unsupported  gate  while  a  hinge-pin  connection 
simulated  the  translational  restraint  of  the  base  hinges  and  prop  connection. 
Analytical  results  shown  in  Table  3  clearly  show  the  increasing  sensitivity  of  the 
higher  modes  to  the  boundary  conditions  of  the  wicket.  An  unchanged  first 
bending  mode,  however,  indicates  its  independence  of  the  boundary  conditions. 
Modeling  of  boundary  conditions  was  thus  recognized  as  one  of  the  most  critical 
tasks  during  the  development  phase  of  this  model  study. 

Sensitivities  of  the  boundary  conditions,  on  the  global  modes  of  the 
wicket,  were  throughly  examined  during  the  experimental  investigation  of  the  1 :5- 
scale  physical  model  (Chowdhury,  Hall,  and  Davis  in  preparation).  Experimental 
results  indicated  that  vibrational  modes  above  1 00  Hz  were  closely  interacted  with 
the  prop  dynamics  and  were  most  sensitive  to  the  prop-rod  connection 
mechanisms.  Figure  7  shows  the  variation  in  the  modal  peaks  for  the  wicket  when 
scanned  from  upstream  and  downstream  sides  of  the  flume.  A  change  in  the 
frequency-response  function  (FRF)  indicates  that  the  system  response  for  the  gate 
was  sensitive  to  the  direction  of  the  driving  force.  Such  a  change  in  the  peak 
resonant  frequency  resulted  primarily  due  to  the  nonlinear  behavior  of  the  "no¬ 
tension"  supporting  mechanisms  of  the  prop  rod  during  the  excitation  of  the 
wicket.  Also  the  uncertain  orientation  of  the  clevis  connection  at  the  top  end  of 
the  prop  rod  introduced  non-linearity  into  the  system  performance.  Experimental 
modal  resonant  frequencies  for  upstream  and  downstream  scanning  are  also 
presented  with  the  analytical  results  in  Table  3. 
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Table  3 

Analytical  and  Experimental  Modal  Frequencies 


Frequency  (Hz) 

Mode 

Analytical 

Experimental 

Description 

Boundary 

Conditions 

Scanning  Laser  Positions 

Free- 

free 

Hinge- 

pin 

Up-stream 
face  scanned 

Down-stream 
face  scanned 

Mode  1 

49.11 

43.11 

Mode  2 
(bending) 

57.74 

57.74 

62.26 

58.00 

Mode  3 
(torsion) 

149.09 

97.99 

95.57 

91.40 

Mode  4 
(bending) 

251.25 

113.91 

171.70 

141.20 

Mode  5 
(torsion) 

314.82 

133.50 

202.01 

232.60 

upstream  scarniing 
dowiLStream  scanning 


|i»-  'V 

Frequency'  (Uz) 
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Analytically  computed  first  bending  and  torsional  modes  occurred  at 
frequencies  of  57.74  Hz  and  97.99  Hz,  respectively,  and  second  bending  and 
torsional  modes  occurred  at  1 13.91  Hz  and  133.50  Hz,  when  the  supports  were 
modeled  as  hinge-pin.  Vibrational  shapes  for  these  modes  are  shown  in  Figure  8. 
These  fi-equencies  are  compared  with  experimental  results  in  Table  3.  A 
comparison  of  the  analytical  and  experimental  modal  frequencies  indicates  that 
the  first  elastic  mode  below  1 00  Hz,  for  the  hinge-pin  condition,  approached  the 
experimental  result  despite  the  fact  that  these  modes  were  not  absolutely  time- 
invariant.  However,  predicted  frequencies  for  the  other  structural  modes  did  not 
quite  match  the  experimental  results. 

As  indicated  in  the  experimental  results,  the  higher  modes  over  1 00  Hz 
would  be  very  difficult  to  match  during  the  modal  updating  due  to  the  intrinsic 
uncertainty  and  non-linearity  of  the  system.  A  number  of  higher  unstable  wicket 
modes  varied  greatly  due  to  their  interaction  with  the  localized  prop  modes. 
Although  the  peak  frequencies  differed,  the  first  five  mode  shapes  were  found 
most  stable  and  repeatable  during  the  experiments.  For  these  modes,  mode  shapes 
did  not  differ,  however,  the  modal  frequencies  changed  during  the  experiments. 
Modal  assurance  criteria  (MAC)  matrix  for  the  first  five  mode  shapes  obtained 
from  different  experiments  showed  a  very  good  agreement  between  the  identical 
modes.  Therefore,  for  the  remaining  of  the  study  only  the  first  five  global  modes 
considered  to  be  most  stable,  repeatable,  and  reliable,  were  used  during  the  modal 
correlation  and  updating  process. 

Modeling  of  the  uncertain  orientation  of  the  clevis  connection  at  the  top 
of  the  prop  rod  is  beyond  the  scope  of  the  present  research.  Such  an  intrinsic 
uncertain  random  behavior  could  be  modeled  using  the  stochastic  FE  method 
(Chang  1993).  The  bottom  end  of  the  prop  rod  in  its  locked  position,  provides 
restraint  only  along  a  direction  away  from  the  gate. 
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Figure  8.  Calculated  mode  shap)es  for  the  initial  FE  model 
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static  response 

A  preliminary  comparison  of  the  static  test  results  with  that  of  the 
analytical  results  was  conducted  to  measure  the  adequacy  of  the  static  stiffness 
used  for  the  analytical  model.  The  static  analysis  modeled  the  test  in  which  a  500- 
Ibf  static  load  was  distributed  across  the  top  surface  of  the  gate  as  shown  in  Figure 
9.  This  loading  condition  corresponded  to  test  load  case  1 .  In  which  a  cylindrical 
disk  mounted  on  a  12- in.  by  4.56-in.  steel  plate  was  placed  on  the  top  edge  of  the 
skin  plate.  Test  load  case  1  applied  a  500-lbf  load  on  a  center-line  axis  passing 
through  the  top  two  gauge  locations.  Five  extensometers  attached  to  the  gate  were 
used  to  record  the  linear  deflection  at  gauge  locations.  Displacement  gauge 
locations  and  FE  model  nodes  at  those  locations  are  also  shown  in  Figure  9.  The 
test  gate  was  supported  by  the  lifting  cylinder  at  -1.1 -deg  down  from  the 
horizontal  axis.  At  its  resting  position,  the  gate  lies  at  -3.3-deg  down  the 
horizontal  axis. 


Two  values  of  spring  stiffness  were  used  for  the  hydraulic  cylinder. 
Analysis  1  used  the  initial  value  of  94,035  Ibf/in.  based  on  the  stiffness  of  the  oil 
column.  A  second  value,  15,1761.6  Ibf/in.,  was  determined  by  calculating  the 
reaction  at  the  ball  support  for  zero  vertical  displacement  and  dividing  this  value 
by  the  test  displacement  of  0.008  in.  at  this  point.  The  second  case  will  be 
identified  as  Analysis  2.  Calculated  reactions  using  the  unyielding  cylinder 
support  for  this  preliminary  analysis  were  -714  Ibf  at  the  gate  hinges  and  1,214  Ibf 
at  the  ball  support.  Note  that  the  stiffness  of  the  mid-point  support,  provided  by 
the  hydraulic  cylinder,  was  time-dependent  and  caused  the  gate  to  deflect 
independent  of  the  loading  situation.  A  time-depended  deformation  of  the 
hydraulic  cylinder  was  possible  due  to  relaxation  of  the  uncertain  amount  of  air- 
pockets  trapped  inside  the  hydraulic  chamber.  As  a  result  of  such  experimental 
shortcomings  it  was  impossible  to  measure  the  true  deflection  of  the  gauge 
locations  during  the  static  test. 


Figure  9.  Location  of  distributed  loading  and  displacement  gagues  for 
static  experiment 
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In  the  analysis  using  the  94,035  Ibf/in.  spring  stiffness,  displacements  at 
nodes  at  the  center  and  base  of  the  gate,  corresponding  to  gauges  3  through  5, 
were  greater  than  test  displacements,  indicating  that  the  spring  stiffness  was  too 
low.  Displacements  at  these  locations  in  the  second  analysis  were  close  to  test 
displacements,  indicating  a  more  nearly  correct  spring  stiffness.  However, 
displacements  at  the  top  of  the  gate  in  the  second  analysis  were  only  about  77 
percent  of  the  test  displacements.  These  differences  could  be  attributed  to  the 
shortcomings  of  the  experimental  results  mentioned  earlier.  Test  displacements 
and  displacements  for  the  two  analyses  are  compared  in  Table  4.  Reactions  at  the 
gate  hinges  for  each  analysis  are  shown  in  Table  5. 


Table  4 

Experimental  and  ABAQUS  Results  for  Static  Test  1 

Gauge  No. 

Node 

Displacements  (in.) 

Test 

Analysis  1 

Analysis  2 

Results 

Results 

1 

10131 

0.048 

0.0448 

0.03707 

2 

9367 

0.048 

0.0448 

0.03705 

3 

5105 

0.008 

0.0112 

0.00797 

4 

3489 

0.003 

0.0049 

0.00313 

5 

3458 

0.003 

0.0048 

0.00308 

Table  5 

Gate  Hinge  Reactions,  Static  Analysis 

Node 

Reactions  (Ibf) 

Test 

Analysis  1 

Analysis  2 

428 

-218.2 

-275.9 

453 

-203.1 

-262.8 
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Revised  FE  Model 


Since  the  FE  model  duplicated  the  1:5  test  gate  model  materials, 
densities,  and  sections,  a  high  degree  of  confidence  could  be  placed  in  the 
accuracy  of  the  model.  Mass  in  the  FE  model  was  verified  by  comparison 
between  the  ABAQUS  and  actual  values.  The  good  correlation  between  the 
ABAQUS  static  test  predictions  and  the  experimental  data  confirmed  that  the 
gate  stiffness  was  reasonably  accurate. 

However,  the  accuracy  of  the  predicted  vibrational  characteristics  of  an 
elastic  FE  model  is  dependent  on  three  factors:  stiffness,  mass,  and  boundary 
conditions.  In  an  attempt  to  more  accurately  model  the  1 : 5-scale  physical 
wicket’s  support  conditions,  a  detailed  modeling  of  the  supports  was  undertaken. 
Additional  modifications  were  made  to  simulate  the  surrounding  water  for 
operating  condition  response  spectrum  analyses  and  to  provide  damping 
coefficients  for  a  time-domain  dynamic  analysis.  The  revised  FE  model  is  shown 
in  Figure  10.  This  figure  shows  a  side  view  of  the  prop-supported  wicket  model 
along  with  the  global  directional  axes.  As  shown  in  the  figure,  the  physical  details 
of  the  wicket  geometry  were  explicitly  discretized  in  the  FE  model.  A  more 
detailed  isometric  view  of  the  FE  model  is  shown  in  Figure  11.  This  figure 
illustrates  the  idealized  boundary  conditions  for  the  base-hinges  and  the  prop 
bottom  end  connection  for  the  wicket  assembly.  Translational  restraints  with 
rotational  flexibility  at  the  mid-point  of  the  hinge  pins  represented  the  actual 
holding  of  the  hinge  pins  within  the  rigid  lugs  mounted  on  the  sill  bearing  frame. 
Translational  restraint  at  the  prop-bottom  end  was  used  to  model  the  fixed-gate 
condition  (65-deg  position)  of  the  prop  blade  while  seated  on  the  hurter  recess  as 
shown  in  Figure  4b. 
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Figure  10.  Revised  FE  model  with  reference  axes  for  tests  and 
analyses 


Modeling  of  supports 

The  gate  base-hinges,  support-rod  hinge  and  the  pillow-block  were 
modeled  using  approximately  300  additional  20-node  brick  and  15-node 
wedge  elements.  Solid  elements  with  no  mass  were  used  to  model  hinge 
openings.  The  prop-  rod  and  hinge  pins  were  modeled  with  2-node  beam 
elements.  Continuity  between  the  3  DOF  solid  elements  and  the  6  DOF 
shell  elements  used  to  model  the  gate  was  imposed  by  the  use  of  multi¬ 
point  constraints  (MFC’s)  at  element  interfaces.  Such  an  MFC  ensures 
proper  moment  transferring  from  the  shell  to  the  interconnected  solid 
elements  (Ridlon  1987).  The  shell-solid  intersection  was  constrained 
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such  that  the  rotational  DOF  at  the  intersecting  nodes  were  compatible. 

The  kinematic  constraints,  MFC  types  SS  LINEAR  and  SS  BILINEAR, 
were  used  depending  upon  the  location  of  nodes  and  the  type  of  elements 
intersecting  each  other.  At  the  prop-rod  hinge,  separate  nodes  were  used 
for  the  prop-pin  beam  elements  and  the  no-mass  solid  elements  at  the 
openings,  and  by  pinned  supports  at  the  hinge-pin  centers.  Figure  12  and 
13  illustrate  the  FE  idealization  of  the  base-hinge  and  prop-end 
cormections,  respectively.  The  prop-rod  congruent  nodes  were  tied  with  pin- 
type  MFC’s  to  allow  rotation  between  the  beam  elements  and  the  solid 
elements.  At  the  gate  hinges,  rotation  was  allowed  but  the  end  was 
restrained  in  the  y-direction  shown  in  Figure  10  and  in  the  x-direction,  along  the 
width  of  the  gate.  For  the  response  spectrum  calculations,  two  conditions  were 
simulated  for  the  restraint  in  the  z-direction;  full  restraint  and  no  restraint.  To 
simulate  the  slotted  connection  used  at  the  prop-rod  end  in  the  1 :5  physical 
model,  only  negative  z-direction  displacements  were  restrained  in  the  dynamic 
time-history  analysis. 


Figure  1 1.  An  illustration  of  boundary  conditions  in  the  FE  model 
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Figure  12.  Analytical  simulation  of  base  hinge-pin  connection 


Figure  13.  Analytical  simulation  of  prop-rod  and  gate  cormection 


2. 
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Added  fluid  masses 


Physical  parameters  obtained  from  the  experimental  results  provided  the 
most  valuable  information  for  updating  the  FE  model  of  the  proposed  wicket. 
These  included  the  static  and  dynamic  parameters  indicative  of  the  actual 
performance  of  the  operating  wicket.  Dynamic  parameters  included  the  dry  and 
wet  gate  modal  parameters  and  the  physical  input  loading  and  output  response 
history  from  the  actual  flow  measurements.  Wet  gate  natural  characteristics 
were  synthesized  to  determine  the  added-mass  contribution  representing  the 
effects  of  water  upon  the  dynamics  of  the  wicket. 

To  account  for  the  effects  of  fluid  on  the  elastic  response  of  the 
structure  under  operating  conditions,  fluid  added  masses  were  distributed  across 
the  top  surface  of  the  FE  model.  Experimental  vibrational  characteristics  for  the 
dry  gate  and  for  the  gate  subjected  to  various  operating  conditions  were  obtained 
from  modal  experiments  using  the  SLDV  system.  Using  the  plate  vibration 
theory,  the  added  masses  for  different  pool  configurations  were  computed  for  the 
operational  wicket  (Chowdhury,  Hall,  and  Davis  in  preparation).  The  added 
mass  (Am)  for  each  mode  was  calculated  from  experimental  modal  frequencies 
and  amplitudes  using  the  formula  shown  below. 


A»i,  =  m 


'  fu 


4ry 


-1 


(1) 


where  f/  is  frequency  in  Hertz  for  mode  i  for  dry  and  wet  conditions  and 
m  is  mass  in  lb-sec2/in.  Experimental  results  showed  that  the  added  mass  varied 
with  mode  and  operating  configuration.  From  an  analytical  point  of  view  a 
single  added-mass,  however,  would  be  convenient  for  computing  global 
response  for  the  wet  wicket.  A  weighted  average  added-mass  was  suggested 
based  on  the  contribution  of  each  modal  peaks  in  the  wet  frequency  response  of 
the  wicket. 


Combined  weighted  mass,  ,  for  modes  1  through  5  was  calculated  using  the 
formula: 


hm 


ZU-Aw,) 


(2) 


where  A/  is  amplitude  at  resonant  frequency  i  for  the  wet  wicket. 

Results  for  the  21-ft  head  condition  are  summarized  in  Table  6.  Note 
that  the  sensitivity  study  of  the  added-mass  conducted  in  the  later  part  of  the 
study  would  determine  the  effects  of  added-mass  upon  the  dynamics  of  the  wet 
gate.  Neglecting  the  contribution  of  second  mode,  in  Table  6,  the  weighted 


Chapter  2  FE  Model  -  Development  Phases 


25 


average  added-mass  could  be  estimated  as  0.4366  Ib-sec^/in. 


Damping 

Damping  has  little  effect  on  predicted  modal  shapes  and  frequencies,  but 
can  affect  results  in  a  dynamic  time-domain  analysis.  In  many  situation, 
however,  damping  has  only  small  effects  on  the  predicted  response,  so  that 
approximate  damping  formulations  are  sufficient  in  the  FE  analysis  (Meyer  and 
Will  1987).  Rayleigh  damping  is  one  of  the  widely  used  models  for  damping 
characteristics  of  an  FE  system.  Rayleigh  damping  coefficients  a  and  P  used  in 
the  dynamic  analysis  were  calculated  from  experimental  data  using  the  formula 
below. 


_ 1 

[2(a{a  +  Pa)f)\  ^ 

where  c/  is  damping  for  mode  i, 

®/  is  27t/T/  and 

T/  is  period  for  mode  i  in  sec. 

Coefficients  for  the  21 -ft  head  operating  condition,  calculated  from 
modes  1  and  3,  were  used  in  the  time-history  analysis.  Experimental  results  for 
the  modal  damping  in  the  wet  case  are  presented  in  Table  6. 


Table  6 

Experimental  Dynamic  Characteristics,  Dry  and  21 -ft  Head  Conditions 
for  the  Prop-Supported  Wicket 

Mode 

fdry  (Hz) 

fwe.  (Hz), 

(percent  damping ) 

A 

^wet 

(in./sAlbf) 

Ani| 

(Ib-sec^/in.) 

1 

43.1 

38.7  (2.6) 

56.3 

0.118168 

2 

58.0 

58.6  (1.9) 

460.4 

0.0 

3 

91.4 

68.6  (3.0) 

477.8 

0.38174 

4 

141.2 

96.8  (3.1) 

147.4 

0.55453 

5 

231.6 

148.47  (2.9) 

203.7 

0.49855 

weighted 

average 

Am„  =  0.3090 
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Verification  of  the  Revised  FE  Model 


Modal  and  dynamic  flow  experiments  provided  experimental  data  for 
verification  of  the  FE  model; 

1 )  The  results  of  modal  testing  under  “dry”  and  operating  conditions 
for  the  prop-supported  wicket  were  compared  with  the  results  of 
response  spectrum  calculations.  Figure  14  shows  a  schematic 
diagram  of  a  typical  experimental  setup  for  the  modal  experiment 
using  the  SLDV  system.  Wet  modal  experimental  data  for  the 
operating  conditions  of  a  21-ft  head  difference,  and  a  three-gate  gap 
(3GG)  flow  configuration  were  used  for  updating  the  FE  model.  A 
3GG  flow  condition  for  the  prop  supported  gate  is  shown  in  Figure 
15. 


Figure  14.  Experimental  modal  setup  schematic  diagram 
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Figure  15.  A  three-gate  gap  flow  configuration 


2)  Dynamic  flow  tests  were  conducted  to  simulate  various  operational 
conditions,  including  varying  heads  of  water  and  gate  positions. 
Dynamic  flow  tests  are  discussed  by  Chowdhury,  Hall,  and  Davis  (in 
preparation).  Although  a  wealth  of  experimental  data  was  available  from 
the  various  tests,  the  estimated  CPU  time  required  for  a  dynamic  FE 
analysis  was  more  than  100  hours,  necessitating  several  restarts.  Due  to 
time  constraints,  a  single  test  was  selected  for  FE  simulation.  Measured 
pressure  distributions  were  used  to  load  the  FE  model,  and  accelerations 
from  the  analysis  were  compared  with  test  accelerometer  records. 


Response  spectrum  results 

Since  variable  restraint  at  the  base  of  the  support  rod  could  not  be 
modeled  in  the  response  spectrum,  two  analyses  for  each  load  condition  were 
run  for  the  dry  and  21 -ft  head  conditions:  full  restraint  in  the  global  z  direction 
and  no  restraint  in  the  global  z  direction.  The  first  test  mode,  rigid-body 
vibration,  could  not  be  reproduced  in  any  analysis  in  which  the  end  of  the 
support  rod  was  fully  restrained.  Analyses  with  no  restraint  produced  rigid-body 
modes  at  somewhat  lower  frequencies  than  the  test  frequencies  and  produced 
second  bending  modes  at  significantly  higher  frequencies  than  the  analyses  with 
restraint.  The  frequencies  of  the  first  bending  and  torsional  modes  and  the 
second  torsional  mode  were  not  affected  by  support  rod  restraint  conditions. 

Predicted  frequencies  for  the  higher  modes  (above  100  Hz)  did  not  agree 
well  with  experimental  results.  For  the  wet  condition,  the  experimental  second 
bending  mode  occurred  at  a  frequency  of  96.8  Hz.  ABAQUS  predicted 
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frequencies  for  the  second  bending  mode  in  the  wet  condition  were  80.035  Hz 
for  restraint  in  3  directions  at  the  end  of  the  support  rod,  and  157.2  Hz  for 
restraint  in  the  global  x  and  y  directions  only.  Two  factors  account  for  the 
inaccuracy  of  the  ABAQUS  predictions.  These  are; 

a.  The  actual  test  condition,  in  which  z  displacement  in  the 
direction  of  water  flow  was  restrained,  lies  somewhere  between 
the  two  modeled  conditions. 

b.  The  added  masses  are  much  larger  for  the  second  bending  and 
torsional  modes  than  for  the  first  3  modes. 


An  additional  analysis  was  run  for  the  wet  condition  without  z  restraint 
at  the  end  of  the  support  rod  and  using  an  added  mass  of  0.5  Ib-sec^/in.  to  reflect 
the  higher  added  mass  of  the  fourth  and  fifth  modes.  This  analysis  produced  a 
second  bending  mode  at  138  Hz.  Obviously  the  addition  of  mass  does  lower  the 
frequency  of  the  second  bending  mode  if  all  other  variables  remain  equal. 
However,  the  inability  to  model  the  actual  support  conditions  in  a  linear  elastic 
analysis  seems  to  be  a  larger  contributor  to  inaccuracies  at  the  higher  modes. 


Frequencies  calculated  for  the  analytical  model  are  presented  in  Tables  7 
and  8.  In  these  tables,  the  sequential  modal  order  for  the  physical  and  analytical 
models  are  shown  in  the  parenthesis.  Analytical  mode  shapes  for  both  restraint 
conditions  of  the  prop-supported  dry  gate  are  shown  in  Figures  16  and  17.  All 
but  the  fourth  experimental  mode  shape  agreed  very  well  with  that  of  the 
analytical  results  for  the  wicket  model  with  no  z-restraint.  Experimental  mode 
shapes  for  the  prop-supported  wicket  are  shown  in  Figure  18.  As  seen  in  this 
figure,  unlike  the  analytical  mode,  the  fourth  experimental  mode  showed  weak 
torsional  characteristics.  The  degree  of  correlation  between  the  corresponding 
anal5tical  and  experimental  mode  shapes  were  computed  using  the  MAC  matrix 
as  discussed  below. 


A  MAC  number  is  an  indication  of  the  degree  of  statistical  correlation 
between  mode  shapes  obtained  from  two  different  sources.  The  MAC  matrix 
provides  a  measure  of  the  least-squares  deviation  of  the  points  of  mode  shape 
from  the  straight  line  correlation  (Allemang  and  Brown  1982)  The  MAC  matrix 
betweeii  corresponding  mode  shapes  obtained  from  source  A,  ^  ^^d  source 
B,  I  [  could  be  defined  by: 


(4) 


Chapter  2  FE  Model  -  Development  Phases 


29 


*  Shaft  supported  wicket  had  a  torsional  mode  at  87.21  Hz 

**  5th  experimental  mode  was  disregarded  because  of  its  inconsistent  behavior 


due  to  interaction  with  prop  dynamics. 
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Figure  16.  Mode  shapes,  dry  condition,  prop-rod  end  restrained  in  x,  y,  and  z 
directions 


Figure  17.  Mode  shapes,  dry  condition,  prop-rod  end  restrained  in  x  and  y 
directions  only 
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Figure  18.  Experimental  dry  mode  shapes  for  the  prop-supported  physical 
model 


The  MAC  matrix  for  the  dry-gate  mode  shapes  obtained  from  the 
analytical  (no  Z  restraint)  and  experimental  sources  is  presented  in  Table  9. 
These  results  indicate  that  a  high  degree  of  correlation  existed  between  the 
modes  compared  in  this  table.  Similarity  of  multiple  modes,  particularly  the 
first  two,  is  evident  from  this  analysis.  Experimental  investigation  indicated  that 
coupled  modes  in  the  frequency  range  100  to  200  Hz  are  time  variant  and  highly 
dependent  upon  the  prop  dynamics.  It  is,  thus,  highly  unlikely  that  the 
correlation  between  the  analytical  and  experimental  results  will  be  enhanced 
without  incorporating  stochastic  FE  analysis  or  nonlinear  system  identification 
techniques  into  the  current  approach.  Moreover,  the  participation  factor 
obtained  from  the  analytical  results  indicated  that  the  most  influential  mode 
(mode  #5)  that  contributed  significantly  to  the  dynamics  of  the  dry  gate  had  a 
quite  high  correlation  between  the  experimental  and  analytical  shapes.  A 
participation  factor  indicates  how  strongly  motion  in  the  global  directions  is 
represented  in  the  eigenvector  of  that  mode  (HKS  1 994).  For  all  practical 
purposes  the  analytical  model  is  considered  satisfactory  to  represent  the  physical 
dry  modal  behavior  of  the  wicket  model. 

The  MAC  relations  between  the  corresponding  wet  experimental  and 
analytical  mode  shapes  are  shown  in  Table  10.  This  table  confirms  the 
similarity  of  the  identical  modes  presented  in  Table  8.  The  physical  model 
showed  a  number  of  deceptive  modes  that  did  not  exist  in  the  analytical  model. 
Mode  numbers  4  and  5  in  the  experimental  model,  for  instance,  were  the  biased 
modes  due  to  nonlinear  interaction  of  the  prop  dynamics.  A  relatively  higher 
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correlation  value  for  the  6th  mode  was  regarded  as  the  distinguishing  criteria  for 
selecting  the  most  probable  similar  mode  that  existed  between  the  analytical  and 
experimental  models. 


Table  9 

Correlation  Matrix  for  the  Prop-Supported  Dry  Wicket  Modes 
Obtained  From  the  Analytical  (No  Z-Restraint  Case)  and  the 
Experimental  Sources _ 


Row  Source:  Laser_down,  Column  Source:  NEW_6512 

Unlike  the  1  st  rotational  experimental  mode,  the  1  st  analytical  mode  had 
a  bouncing  pattern  as  indicated  in  the  MAC  relation  table.  Such  a  change  in  the 
1  st  rigid  body  mode  from  bouncing  to  rotation  could  be  attributed  to  the 
boundary  effects  on  the  physical  model.  Experimental  results  indicated  that  the 
actual  pattern  of  the  1st  mode  depended  primarily  upon  the  orientation  of  the 
prop-rod  bearing  housing  connection  (top-end)  with  respect  to  the  prop- 
connection  pin  that  was  mounted  tangential  to  the  gate  transverse  axis.  A  slight 
tilting  of  the  bearing  housing,  which  is  quite  obvious  during  the  wet  experiment, 
could  definitely  change  the  resistance  mechanisms  of  the  prop-supported  gate. 
Experimental  modal  results  for  the  wet  and  dry  wicket  indicated  a  bouncing  1st 
mode  for  a  shaft  supported  wicket  rather  than  the  rotational  pattern  obtained  for 


Analytical 
( no  Z  restraint) 

Modal  Order 
(Row) 

1 

(30.2  Hz) 

2 

(57.7  Hz) 

3 

(93.9  Hz) 

4 

(128.0  Hz) 

5 

(205.9  Hz) 

6 

(212.8  Hz) 
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the  prop-supported  gate.  A  one-gate  gap  wet  condition  also  showed  1  st 
bouncing  mode  instead  of  rotational  one  for  the  prop-supported  wicket.  These 
results  support  the  above  conclusion  that  the  interchangeability  of  the  1  st  mode 
is  dependent  upon  the  instantaneous  boundary  condition  of  the  prop  top-end 
connection  and  the  uneven  water  mass  distribution,  if  any,  on  the  submerged 
surface  of  the  wicket.  Realistically,  however,  a  nonuniform  added  mass 
distribution  from  a  uniform  flow  pattern  is  nonobvious.  Considering  such  a 
strenuous  boundary  incompatibility  between  the  physical  and  analytical  models, 
reasonably  good  MAC  numbers  between  the  corresponding  most  contributory 
modes  indicate  a  good  analytical  representation  of  the  wet  physical  model.  As 
discussed  below,  the  strong  modes  that  are  most  contributory  to  the  dynamics  of 
the  wet  wicket  are  identified  from  the  modal  participation  factor  obtained  from 
the  FE  analysis.  In  the  wet  analysis,  the  5th  mode  is  the  most  active  mode  that 
has  the  largest  modal  participation  in  excitation  along  the  z-direction  (normal  to 
the  gate  surface). 


Table  10 

Correlation  Matrix  for  the  Prop-Supported  Wet  Wicket  Modes 
Obtained  From  the  Analytical  (No-Z  Restraint  Case)  and  the 
Experimental  Sources 

Analytical 
( no  Z  restraint] 

Modal  Order 
(Row) 

Experimental  Modal  Order  (Column)  for  21 -ft,  3GG  configuration 

1 

(38.9  Hz) 

2 

(58.6  Hz) 

3 

(68.6  Hz) 

4 

(96.8  Hz) 

5 

(148.4Hz) 

6 

(163.2Hz) 

1 

(22.9  Hz) 

0.890 

2 

(57.7  Hz) 

0.59 

3 

(72.7  Hz) 

0.930 

0.981 

0.506 

4 

(83.3  Hz) 

0.924 

0.979 

0.533 

5 

(157.2  Hz) 

0.826 

6 

(200.0  Hz) 

0.477 

0.699 

Row  Source:  Laser_21_ft_2_g 

•a 

3 

O 

o 

umn  Source:  NEW_6512w 
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Modal  participation  factors  for  the  wet  analyses  show  that  in  the  analysis 
with  z  restraint,  mode  2,  the  first  torsional  mode,  is  the  largest  contributor  to 
gate  behavior.  Mode  6,  a  third  bending  mode  that  occurs  at  218.77  Hz,  is  also  a 
large  contributor.  In  the  analysis  without  z  restraint,  the  first  torsional  mode, 
mode  3,  is  again  the  largest  contributor  to  gate  behavior,  but  the  third  bending 
mode  disappears.  In  both  analyses,  the  first  bending  modes  contribute  little  to 
gate  behavior.  Participation  factors  for  the  first  5  modes  for  the  wet  condition 
analyses  are  given  in  Tables  1 1  and  12. 

The  analytically  computed  modal  effective  mass  also  indicated  that  the 
first  five  modes  could  adequately  represent  the  actual  dynamics  of  the  dry  and 
wet  wicket  assembly.  This  conclusion  was  based  on  the  fact  that  the  sum  of  the 
modal  effective  masses  for  the  first  five  modes  along  the  z-kinematic  directions 
was  about  70  percent  of  the  total  mass.  This  suggests  that  the  first  five  modes 
have  significant  participation  in  excitation  along  the  z-direction,  normal  to  the 
gate  surface.  The  higher  modes  contribution  (modes  6  through  10)  was  very 
insignificant,  in  the  order  of  10  which  could  easily  be  neglected  for 
determining  the  total  response  of  the  wicket  system.  The  sum  total  of  modal 
effective  masses  and  the  contribution  along  the  kinematically  restrained  degrees 
of  freedom  should  provide  the  total  mass  of  the  model  (HKS  1994).  In  this  case, 
30  percent  of  total  mass  remained  fixed  along  the  z-direction  due  to  the 
kinematic  restraints  at  the  supports. 


Time-domain  dynamic  response 

Response  spectrum  analyses  indicated  that  gate  behavior  was  well 
modeled  for  the  first  bending  and  torsional  modes,  but  higher  modes  were  more 
directly  influenced  by  prop-rod  end  conditions  and  were  less  well  modeled. 
Based  on  this,  an  initial  decision  was  made  to  attempt  to  model  only  frequencies 
less  than  100  Hz  in  the  dynamic  analysis.  The  maximum  time  step  that  can  be 
expected  to  allow  excitation  of  the  structure  at  a  frequency  of  1 00  Hz  is 
approximately  2.5  times  the  period,  or  0.004  sec.  Analysis  of  the  frequency 
content  of  the  acceleration  records  for  the  various  experiments  revealed  that  for  a 
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given  experiment,  the  frequency  content  of  each  third  of  the  data  was  similar,  and 
only  one-third  of  the  test  records  were  used  for  loading  the  model  and  evaluating 
resulting  accelerations.  During  the  experiment  the  accelerometer  data  were 
filtered  at  500  Hz,  thus  the  resampled  data  used  in  the  analysis,  having  a  third  of 


Table  11 

Modal  Participation  Factors,  “Wet”  Analysis,  Z  Restraint 

Mode 

Participation  factors,  analysis  with  z  restraint  on  prop  rod 

x-comp 

y-comp 

z-comp 

x-rot 

y-rot 

1 

0.0005 

0.01723 

0.04371 

2.8714 

-0.0085 

2 

438.79 

3.2717 

-0.48626 

98.602 

-8005.6 

3 

-0.002 

0.56161 

1.4385 

58.422 

0.11666 

4 

245.55 

0.32975 

-6.3558 

-168.93 

-1675.4 

5 

0.22362 

0.01294 

-0.02699 

-0.06653 

-5.6198 

Table  12 

Modal  Participation  Factors,  “Wet”  Analysis,  no  Z  Restraint 

Mode 

Participation  factors,  analysis  with  no  z  restraint  on  prop  rod 

x-comp 

y-comp 

z-comp 

x-rot 

y-rot 

1 

0.0004 

0.61672 

1.3814 

43.895 

0.07844 

2 

0.0008 

-0.01688 

-0.03679 

0.59786 

-0.01342 

3 

438.95 

2.9482 

-1.29 

67.969 

-8008.5 

4 

244.57 

2.097 

-1.7118 

25.925 

-1668.4 

5 

2.0897 

54.913 

116.78 

-945.3 

-22.216 
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the  total  frequency  range  of  the  experiments,  almost  contained  the  first  five  wet 
modes  of  the  physical  model.  Test  accelerometer  locations  and  corresponding 
nodes  are  shown  in  Figure  19.  In  this  figure,  the  accelerometer  locations  are 
designated  by  the  suffix  A/,  where  /  =  1  to  9. 
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Figure  19.  Location  of  FE  nodes  and  accelerometers 


The  dynamic  analysis  modeled  experiment  number  163,  with  the  gate  in 
the  65-deg  position  and  a  21-ft  head  of  water.  A  static  loading  equal  to  an 
average  of  the  test  pressure-time  histories  for  the  first  few  load  increments  was 
used  as  in  initial  loading  for  the  gate.  Displacement  in  the  direction  opposite  the 
water  flow  at  the  base  of  the  support  rod  was  allowed  by  providing  a  gap  element 
between  the  node  at  the  end  of  the  support  rod  and  a  fixed  node  at  the  same 
location.  This  allowed  free  movement  of  the  support  rod  end  away  from  the  fixed 
node,  but  prevented  further  movement  in  the  direction  of  flow  once  the  gap  was 
closed.  Test  downstream  pressures  were  subtracted  from  upstream  pressures  to 
give  a  net  pressure  at  each  gauge  location,  and  pressures  for  loading  the  plate 
surface  were  calculated  by  interpolating  between  test  data  for  the  9  gauge 
locations,  shown  in  Figure  20.  Thus,  the  net  pressure  envelope  comprising  the 
nine  gauge  readings  provided  the  flow-induced  pressure  loads  for  the  analytical 
model. 
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Figure  20.  Pressure  gauge  locations  and  column  lines  for  displaying  pressure 
readings 


An  initial  analysis  using  approximately  0.5  sec  of  the  pressure  history 
showed  that  peak  accelerations  were  calculated  at  each  0.004  sec  increment,  and 
that  times  for  these  peaks  did  not  exactly  correspond  to  those  in  the  test 
acceleration  records.  This  indicated  that  the  0.004-sec  time  step  was  not  small 
enough  to  capture  all  of  the  response,  and  the  analysis  was  reinitiated  with  a 
0.002-sec  time  step.  The  analysis  was  then  conducted  for  5  sec,  requiring 
approximately  120  hr  of  CPU  time  on  the  CRAY  C-90.  Calculated  and 
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experimental  acceleration  records  are  shown  for  eight  locations  in  Figures  21 
through  28.  Calculated  acceleration  records  were  not  saved  at  node  745 1 . 

At  all  of  these  locations,  the  abrupt  application  of  dynamic  loads  to  a 
model  with  no  initial  accelerations  resulted  in  large  initial  accelerations  that  were 
damped  to  more  consistent  values  by  approximately  0.05  to  0.15  sec.  Because  of 
this,  the  first  0. 1 5  sec  of  data  were  deleted  prior  to  evaluating  the  modal  content  of 
the  ABAQUS  acceleration  records.  The  resulting  records  shifted  to  an  initial 
value  of  zero  are  the  records  shown  in  Figures  21  through  28.  Data  at  nodes 
7451,  corresponding  to  gauge  5,  was  not  saved. 

In  the  acceleration  records  for  nodes  near  the  top  of  the  gate,  shown  in 
Figures  21  through  23,  and  at  nodes  near  the  bottom  of  the  gate,  shown  in  Figures 
26  through  28,  the  predicted  magnitudes  of  negative  accelerations  were  close  to 
the  measured  accelerations,  while  positive  accelerations  were  under  predicted.  At 
nodes  near  the  center  of  the  gate,  shown  in  Figures  24  and  25,  negative 
accelerations  were  under  predicted. 

Preliminary  fast  Fourier  transforms  (FFT)  of  the  calculated  and  test 
accelerations  at  the  eight  locations  are  shown  in  Figures  29,  3 1,  and  32.  Peak 
amplitudes  at  60  Hz,  120  Hz,  and  240  Hz,  seen  in  Figures  29  and  31,  appear  in  all 
of  the  test  record  FFTs  and  are  a  function  of  variations  in  electrical  frequency. 
These  frequencies  were  embedded  in  the  pressure  records  used  to  load  the 
analytical  model  and  also  appear  in  the  ABAQUS  records.  Since  amplitudes  were 
low  except  at  these  peaks,  only  amplitudes  through  0.5  have  been  plotted  on  the 
ordinate  in  FFT  plots  shown  in  Figures  3 1  and  32.  Only  frequencies  through  150 
Hz  are  plotted,  since  the  structural  frequencies  of  interest  lie  in  this  range. 
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Z  Acceleration  (in./sec^) 


Figure  28.  Z  accelerations,  node  2769  and  gauge  9 
Chapter  2  FE  Model  -  Development  Phases 


The  frequency  content  of  test  and  calculated  records  is  also  affected  by 
the  dynamic  behavior  of  the  prop  rod.  The  frequency  response  function  for  the 
prop  rod  in  the  21-ft  head  condition  is  shown  in  Figure  30.  Distinct  prop-rod 
modes  occur  at  38,  58,  70,  and  115  Hz.  Modes  at  higher  frequencies  are  not 
distinct  and  are  probably  due  to  interaction  between  the  prop  rod  and  the  gate, 
further  complicating  the  problem  of  correctly  modeling  prop-rod  boundary 
conditions.  In  the  response  spectrum  results,  the  frequencies  for  the  second 
torsional  mode  in  the  analyses  with  no  Z  restraint  dropped  below  those  of  the 
second  bending  mode.  This  shows  the  heavy  dependence  of  gate  dynamic 
response  on  prop-rod  response  in  the  higher  modes. 


Figure  29.  FFT  of  ABAQUS  and  test  records, 
gauge  8  and  node  2858 


The  FFTs  of  the  test  z  acceleration  records  are  shown  in  Figure  3 1 .  Peaks 
near  47  Hz  can  be  seen  in  all  records  except  those  at  the  gate  centerline.  Peaks 
between  52  and  56  can  be  observed  in  all  records  except  for  those  for  gauges  7 
through  9  at  the  bottom  of  the  gate.  Peaks  between  63  and  73  Hz,  at 
approximately  148  Hz,  and  at  approximately  163  Hz  are  evident  in  all  gauges 
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except  those  at  the  gate  centerline,  gauges  2,  5,  and  8.  Peaks  at  approximately  99 
Hz  and  195  Hz  are  evident  in  all  gauge  records.  Any  structural  response  near 
120  Hz  is  obscured  by  peaks  due  to  electrical  fluctuations.  The  first  4 
frequencies,  47  Hz,  52  to  56  Hz,  63  to  73  Hz,  and  99  Hz  may  correspond  to  the 
first  4  mode  shapes,  although  the  lack  of  a  47-Hz  peak  at  centerline  gauges  would 
seem  to  indicate  a  torsional  motion  at  this  frequency.  Peaks  at  higher  frequencies 
are  less  defined,  and  except  for  the  peaks  at  163  Hz,  do  not  correspond  to  any 
distinct  mode  shape  and  may  be  due  to  gate-prop  rod  interaction. 


Frequency  response  function 


Figure  30.  Frequency  response  function  for  prop  rod 


The  frequency  content  of  ABAQUS  acceleration  records  for  nodes 
located  near  8  of  the  9  gauges  is  shown  in  Figure  32.  Accelerations  at  node 
7451,  near  the  center  of  the  gate,  were  not  saved.  Weak  peaks  near  47  and  55  Hz 
and  strong  peaks  near  73  Hz  are  visible  in  records  for  nodes  near  the  top  and 
center  of  the  gate  and  may  correspond  to  the  first  3  mode  shapes.  However, 
predicted  amplitudes  for  the  first  torsional  mode,  at  approximately  73  Hz,  are 
much  larger  than  test  amplitudes.  Peaks  near  99  Hz  are  evident  for  the  3  records 
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at  the  base  of  the  gate  only.  Peak  amplitudes  also  are  found  near  90  Hz,  105  Hz, 
134  Hz,  150  Hz,  and  159  Hz.  These  frequencies  do  not  correspond  to  any  mode 
defined  by  the  response  spectrum  analyses.  The  appearance  these  indistinct  higher 
modes  and  the  low  amplitude  of  the  rigid-body  response,  which  was  heavily 
affected  by  prop-rod  boundary  conditions  in  the  response  spectrum  results,  can  be 
attributed  to  prop-rod  response  and  interaction  between  the  gate  and  the  prop  rod. 

Plots  of  test  amplitude  versus  calculated  amplitude  for  frequencies 
ranging  from  0  to  150  Hz  are  shown  in  Figure  33.  Amplitudes  of  frequencies 
between  115  and  130  were  zeroed  for  this  plot,  since  amplitudes  in  this  range  of 
frequencies  were  dominated  by  the  effects  of  electrical  fluctuation  and  their 
inclusion  resulted  in  much  wider  scatter.  Correlations  are  better  for  nodes  near 
the  top  and  bottom  of  the  gate,  and  poor  for  nodes  located  at  the  center  row  of 
gauges. 


Since  analysis  predictions  were  heavily  dependent  on  prop-rod  support 
conditions,  an  additional  analysis  was  performed.  In  the  second  analysis,  a 
minimum  gap  of  0.01  inch  was  maintained  between  the  prop-rod  end  and  the 
fixed  node.  The  purpose  of  this  analysis  was  to  provide  a  less  stiff  support  at  the 
prop-rod  end.  This  analysis  was  only  continued  for  2.5  sec.  Plots  of  z 
accelerations  at  3  locations  are  shown  in  Figures  34  through  36.  As  in  the 
previous  analysis,  predicted  negative  accelerations  were  similar  to  test  negative 
accelerations  at  the  top  and  bottom  of  the  gate,  and  predicted  positive 
accelerations  were  low.  However,  the  frequency  content  of  the  acceleration 
records,  shown  in  Figure  37,  was  somewhat  different  from  that  of  the  previous 
analysis.  Additional  peaks  between  20  and  30  Hz  and  at  frequencies  higher  than 
80  Hz  confirm  that  gate-prop  rod  interaction  affects  predictions  in  these 
frequency  ranges. 

Test  reactions  were  only  available  at  the  gate  hinges  due  to  failure  of  the 
load  cells  at  the  prop-rod  end.  Global  y  and  z  reactions  at  the  hinges  from  the  two 
ABAQUS  analyses  are  compared  with  test  reactions  in  Figure  38,  ABAQUS 
reactions  at  the  prop-rod  end  are  shown  in  Figure  39,  and  total  global  y  and  z 
reactions  from  the  two  analyses  are  shown  in  Figure  40.  For  purposes  of 
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comparison,  gate  experimental  y  and  z  forces,  shown  in  Figure  10,  have  been 
converted  to  the  ABAQUS  global  axes.  Analysis  boundary  conditions  are  shown 
in  Figure  41. 

In  both  analyses  the  gate  was  thrown  forward  by  the  water  loading, 
causing  the  majority  of  the  load  to  be  carried  by  the  support  rod.  In  the  first 
analysis,  designated  ABl  in  Figures  38  through  41,  some  vertical  load  was 
distributed  to  the  gauge  hinges,  resulting  in  small  positive  vertical  reactions  at  the 
hinges.  In  the  second  analysis,  AB2,  uplift  occurred  at  the  gauge  hinges. 

However,  total  vertical  and  horizontal  reactions  for  the  two  analyses  were  the 
same,  as  can  be  seen  in  Figure  40.  Since  reactions  at  the  gate  hinges  are  so 
sensitive  to  boundary  conditions  at  the  end  of  the  support  rod,  it  may  not  be 
practical  to  attempt  to  model  test  reactions  exactly.  However,  the  two  ABAQUS 
analyses  bounded  the  measured  reactions  and  produced  adequate  values  for  design 
of  the  supports. 

Maximum  principal  stresses  were  low  throughout  the  analyses,  with  the 
largest  stresses  occurring  at  the  gate  hinges  and  the  support  rod  hinge.  Typical 
maximum  principal  stresses  on  the  top  and  bottom  surface  of  the  top  plate  are 
shown  in  Figure  42.  Maximum  principal  stresses  on  the  outer  surface  of  the  gate 
hinges  are  shown  in  Figure  43. 
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(Amplitudes  are  in  arbitrary  units) 
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(Amplitudes  are  in  arbitrary  units) 


Figure  35.  Z  acceleration,  node  7188  and  gauge  4,  analysis  2 
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Amplitude 


Figure  36.  Z  acceleration,  node  2813  and  gauge  7,  analysis  2 


Figure  37.  Frequency  content  of  Z  acceleration  records  for  3  nodes,  analysis  2 
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Figure  38.  Measured  and  calculated  hinge  reactions 


Figure  39.  Calculated  prop-rod  reactions. 


Figure  40.  ABAQUS  total  reactions,  ABl  and  AB2 
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3  FE  Model-  Applications 


Background 

As  discussed  in  the  previous  chapter,  predicted  response  from  an 
updated  FE  model  is  still  an  approximation  of  the  actual  response.  Uncertain 
loading  and  boundary  conditions  make  it  nearly  impossible  to  simulate  the  true 
response  using  a  deterministic  FE  model.  For  all  practical  purposes,  however, 
an  updated  analytical  FE  model  is  a  flexible  tool  which  could  be  used  reliably 
and  confidently  to  evaluate  the  performance  of  prototype  structure.  One  of  the 
most  important  advantages  of  such  an  FE  model  is  in  the  identification  of 
sensitive  parameters  which  could  lead  to  determine  the  bounds  of  structural 
responses.  Bounds  of  structural  responses  could  effectively  indicate  the 
adequacy  of  the  structure  for  extreme  operating  and  sensitive  boundary 
conditions. 

In  this  chapter  the  applications  of  the  FE  model  in  determining  sensor 
location  using  pre-experiment  FE  analysis,  and  the  results  of  the  fatigue  study  to 
determine  the  effects  of  cyclic  loading,  and  the  upper  bounds  of  wicket  response 
due  to  most  probable  limits  of  boundary  and  loading  conditions  are  discussed. 


Sensitivity  Study 

Static  experimental  results  were  used  to  study  the  sensitivity  of  different 
physical  parameters  and  boundary  conditions  of  the  wicket  assembly.  This 
paramatric  study  was  conducted  using  the  preliminary  FE  model  to  identify  the 
most  sensible  parameters  of  the  wicket  geometry.  As  discussed  previously,  no 
supporting  devices  such  as  base  hinges  and  prop  connections  were  modeled  or 
included  in  the  preliminary  FE  model.  In  this  model,  nodes  located  at  the 
supports  were  restrained  to  simulate  the  boundary  conditions  of  the  actual 
structure.  Hinge  bases  were  restrained  for  translational  motions  by  using  pin 
connections  in  the  FE  model.  The  cylinder  connection  was  modeled  in  two 
different  ways  to  measure  the  sensitivity  of  predictions  to  the  modeling  of  the 
supporting  devices.  Two  different  types  of  boundary  conditions  were  considered 
in  the  modeling  technique.  These  included  a  roller  (Z-translational  restraint)  and 
a  spring  connection.  A  roller  connection  represented  the  actual  holding 
mechanism  of  the  ball-cup  support  in  which  the  mounting  ball  connected  to  the 
gate  frame  had  rested  on  the  cup  connected  to  the  shaft  of  the  hoisting  cylinder. 

The  spring  model  simulated  the  spring  action  of  the  hydraulic  device  which  had 
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been  deformed  due  to  compressibility  of  the  air-entrapped  inside  the  fluid 
chamber. 

In  the  preliminary  FE  model,  the  span  length  between  the  base  hinge  and 
ball-cup  support  was  slightly  shorter  than  in  the  physical  model.  This  difference 
existed  because  the  base  hinges,  which  were  2.6  in.  away  from  the  base  of  the 
gate,  were  not  modeled.  Note  that  in  the  preliminary  model,  restrained  nodes 
were  attached  to  the  base  of  the  hinge  mounting  plate  on  the  extreme  end  of  the 
gate  (Figure  6).  In  the  FE  model,  the  roller  or  spring  connection  used  for 
representing  the  cylinder  support  had  been  set  on  nodes  near  the  ball-cup 
connection.  Considering  the  availability  of  nodes,  two  distances  of  21.4  in.  and 
22.3  in.  were  found  for  simulating  the  proper  span  length  of  the  cup-ball  support 
from  the  base  hinge.  In  the  experimental  model,  the  cylinder's  connection  was 
24  in.  away  from  the  hinge  axis.  Predicted  results  from  the  preliminary  FE 
model,  thus,  required  an  adjustment  to  account  for  the  discrepancy  in  the  span 
length.  An  analytical  adjustment  of  the  FE  results  were  made  by  using  a  ratio  of 
the  beam  deflections  for  actual  (24  in.)  and  simulated  (21 .4  in.  or  22.3  in.)  span 
length.  Using  beam  theory,  correction  factors  for  each  gauge  were  computed 
such  that  the  compensation  for  not  using  correct  span  length  was  included  in  the 
adjusted  results. 

Static  experimental  results  for  three  loading  cases  are  presented  in  Table  13 
below.  These  three  loading  cases  could  be  described  as: 

Load  Case  1  -  in  which  a  500-lbf  load  was  applied  in  between  gauge  #1  and  #2 
(Figure  9  for  gauge  locations). 

Load  Case  2  -  same  500-lbf  load  was  applied  on  the  left  side  of  the  gate, 
directly  over  gauge  #1 

Load  Case  3  -  same  500-lbf  load  was  applied  on  the  right  side  of  the  gate, 
directly  over  gauge  #2 


Table  13 

Measured  Deflection  in  inches 

Load  Case 

Input  Load,  Ibf 

Gauge  #1 

Gauge  #2 

Gauge  #3 

Gauge  #4 

Gauge  #5 

I 

500 

0.048 

0.048 

0.008 

0.003 

0.003 

II 

500 

0,0805 

0.0185 

0.0065 

0.007 

-0.001 

III 

500 

0,0215 

0.0805 

0.0125 

-0.001 

0.00775 

The  500-lbf  load  was  a  cylindrical  disk  mounted  on  a  4.56-in.by  12-in. 
steel  plate  which  provided  a  distributed  load  over  the  contact  surface  with  the 
gate.  In  the  FE  model,  this  load  could  be  regarded  either  as  concentrated  loads 
applied  along  the  center-line  of  contact  or  as  uniformly  distributed  (UDL)  over 
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the  surface  of  contact.  Both  of  these  situations  will  be  examined  in  this  report. 
In  the  first  case,  the  concentrated  load  was  applied  along  nodal  lines 
corresponding  to  the  width  of  the  supporting  plate.  In  the  second  case,  the 
distributed  load  was  applied  over  the  elements  holding  the  supporting  plate. 

Spring  stiffnesses  for  the  lifting  cylinder  were  computed  from  the  static 
deflection  and  are  presented  in  Table  14. 


Table  14 

Spring  Stiffness  for  Boundary  Condition  (Spring)  at  Y=21.4  in. 

Load  Case 

Reaction,  Ibf 

Displacement,  in. 

Spring  Stiffness,  Ibf/in. 

Cone.  Load 

1310 

0.008 

163750 

Dist.  Load 

1285 

0.0065 

197692.31 

Physical  parameters 

As  shown  in  Table  1,  uncertainties  of  material  properties  and  physical 
dimensions  are  the  reality  in  a  physical  model.  To  determine  the  effects  of  these 
variations,  a  series  of  FE  analyses  was  undertaken.  Using  the  preliminary 
wicket  gate  model  and  the  two  load  distributions,  bounding  analyses  were 
conducted  for  four  variables:  span  length,  modulus  of  elasticity,  section 
thickness,  and  prop-support  conditions. 

In  the  first  pair  of  analyses  two  span  lengths  were  considered.  These 
are  22.3  in.  and  21.4  in.  apart,  respectively,  from  the  simulated  base  hinge  to 
the  prop  support.  The  next  two  analysis  modeled  two  values  of  modulus  of 
elasticity.  These  two  values  of  modulus  of  elasticity,  (32,213csi,  and  34,287 
ksi,  respectively),  were  used  based  on  the  experimental  results  presented  in 
Table  1 .  The  first  set  was  based  on  the  minimum  value  of  the  modulus  of 
elasticity  for  the  major  structural  components  of  the  wicket.  And  the  second  set 
was  based  on  the  average  magnitudes  of  the  modulus  of  elasticity  of  the  double 
I-beams.  The  I-beams  were  the  main  structural  component  of  the  wicket  that 
ran  vertically  along  the  length  of  the  gate. 

Different  values  of  thickness  for  the  major  structural  components  were 
also  considered  for  evaluation.  The  minimum  and  average  thickness  for  the 
beams  and  transverse  stiffeners  were  used  in  two  separate  FE  analyses  of  the 
wicket.  Table  15  shows  the  variation  of  thickness  for  the  major  structural 
components.  These  changes  were  considered  to  determine  the  effects  of  the 
uncertainty  of  the  thickness  on  the  structural  performance. 

The  final  pair  of  analysis  simulated  the  boundary  condition  at  the  ball 
support  in  two  ways:  as  a  spring,  and  as  a  roller. 
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Table  15 

Thickness  Variations  in  the  Major  Gate  Components 

Component 

Average  Thickness,  in. 

Minimum  Thickness,  in. 

Beam  Flange 

0.3192 

0.3100 

mMSSsm 

0.1793 

0.1500 

Long  Stiffeners 

0.0750 

0.0700 

The  relative  importance  of  variations  in  each  of  thesepairs  of  analyses 
can  be  seen  in  plots  of  deflections  from  analysis  1  versus  deflection  from 
analysis  2,  shown  in  Figures  44  through  47.  Figure  44  shows  the  effects  of  span 
length  on  the  static  response  of  the  wicket  due  to  two  different  load-modeling 
techniques.  Figure  45  presents  the  relationships  of  gauge  deformations  due  to 
changes  in  the  modulus  of  elasticity.  In  this  figure,  results  corresponding  to  the 
two  different  moduli  are  compared  indicating  the  contribution  of  elasticity  on  the 
overall  stiffness  of  the  wicket.  These  figures  indicate  that  the  stiffness  changes 
due  to  variation  in  span  length  and  material  properties  are  marginal  and 
proportional  when  applied  load  was  treated  as  line-load.  For  the  distributed  case, 
asymmetric  loading  (load  case  2  ),  showed  non-proportionality  in  one  of  the 
gauge  deformation  due  to  changes  in  span  length  and  modulus  of  elasticity. 


Figure  44.  Variation  of  deflection  at  gauge  locations  due  to  changes  in  span  length 
(E  =  32,213  ksi  and  roller  model) 
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Figure  45.  Variation  of  deflection  at  gauge  locations  due  to  changes  in  modulus  of 
elasticity  (  Y  =  21.4  in.  and  roller  model) 


The  effects  of  wicket  component  thickness  on  the  overall  stiffness  of  the 
wicket  are  analyzed  in  Figure  46.  In  this  figure,  deformation  results  for  two 
different  material  thicknesses  are  compared.  As  in  the  previous  cases,  this 
figure  indicates  that  the  ratio  of  deformations  for  different  material  thickness  is 
constant.  Boundary  condition,  however,  appeared  to  be  the  most  sensitive 
parameters  among  the  variables  studied  in  this  investigation.  A  change  in 
boimdary  condition  from  roller  to  spring  showed  a  major  flucmation  in  the  ratio 
of  deformations  as  shown  in  Figure  47.  An  effective  modeling  of  the  boundary 
conditions,  thus,  appeared  to  be  the  most  critical  factor  in  the  FE  analysis  of  the 
Olmsted  wicket.  Nonlinearity  of  prop  support  further  complicated  the 
simulation  process.  As  discussed  in  the  previous  chapter,  “no  tension” 
mechanism  of  the  prop  blade  support  and  random  disturbances  of  the  prop  top 
clevis  introduced  nonlinearity  into  the  structural  system.  This  required  that  the 
FE  model  be  evaluated  for  idealized  boundary  conditions  representing  the  actual 
boundary  limits  of  the  wicket  system  such  that  the  response  limits  could  be 
determined. 
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Figure  46.  Variation  of  deflection  at  gauge  locations  due  to  changes  in  component 
thickness  (E  =  32,213  ksi,  Y  =  21.4  in.,  and  spring  model) 


Figure  47.  Variation  of  deflection  at  gauge  locations  due  to  changes  in  support 
conditions  (E  =  32,213  ksi  and  Y  =  21.4  in.) 

These  results  also  indicate  that  the  best  representation  of  the  static  tests 
in  the  FE  model  is  given  by  the  case  in  which  the  lower  modulus  of  elasticity, 
and  spring  support  were  used.  Adjusted  static  results  using  preliminary  FE 
model  for  load  case  1  are  presented  in  Table  16. 
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Table  16 

Adjusted  Results  for  Modified  Thicknesses  and 

Boundary  Condition  (Spring)  set  at  Y=21.4  in. 

E  =  32,213  E3  psi 

Cone.  Load 

Dist.  Load 

Gauge  No. 

Node  FEM 

th_spr_slowcload1 

th_spr_slowdload1 

1 

7630 

0.04214 

0.062873 

2 

8294 

0.04213 

0.013359 

3 

3408 

0.008188 

0.06688 

4 

2120 

0.004330 

0.006862 

5 

2156 

0.004320 

-0.000354 

Pre-experiment  Investigation 


A  preliminary  FE  analysis  was  used  to  determine  the  best  built-in 
accelerometer  locations  for  dynamic  response  investigation  of  the  wicket  gate. 
The  best  sensor  location  was  defined  as  that  point  at  which  excitation  in  a 
specific  vector  direction  resulted  in  the  best  response  of  all  modes  of  interest. 
This  analysis,  thus,  selected  the  sensor  location  so  that  the  accelerometers  were 
not  placed  on  any  nodal  lines  and  the  resonant  peaks  of  the  FRF  responses  were 
dominant  in  the  synthesized  plots.  Synthesized  FRF  plots  were  generated  from 
the  FE  results  using  the  IDEAS-Masterseries  post-processing  module  (SDRC 
1993).  A  detailed  description  of  the  procedure  and  the  sensor  location  for  the 
experimental  wicket  are  presented  by  Chowdhury,  Hall,  and  Davis  (in 
preparation). 


Fatigue  Analysis 

A  post-experiment  fatigue  investigation  of  the  critical  components  was 
performed  separately  using  the  FE  analysis  of  the  wicket  assembly  (Chowdhury, 
Hall,  and  Davis  in  preparation).  Base  hinges  and  the  shear  pins  shown  in  Figure 
48  are  the  most  critical  components  of  the  wicket  system  that  could  be  subjected 
to  fatigue  failure  due  to  reversed  loading.  A  fatigue  analysis  of  the  most 
sensitive  components  was  conducted  using  the  modified  Goodman  line  method 
which  required  the  computation  of  the  mean  and  alternating  normal  stresses  on 
the  structure  based  on  the  maximum  and  minimum  stress  applied  to  the 
component.  FE  analysis  of  the  components  provided  the  maximum  and 
minimum  stresses  for  the  critical  operating  condition.  The  modified  Goodman 
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line  method  determines  the  usefulness  against  possible  fatigue  failure  of  a 
component  subject  to  a  combination  of  static  and  reversed  dynamic  loading 
(Lindeburg  1994). 


For  the  FE  analysis,  the  critical  normal  loads  were  obtained  from  the 
experiment  conducted  for  the  21 -ft  head  difference,  three-gate  gap  configuration, 
when  the  physical  experimental  gate  was  fixed  at  24  deg  position.  The  predicted 
loads  were  applied  on  a  finite-element  model  of  the  hinge  base  and  the 
maximum  and  minimum  Von-Mises  stresses  were  used  to  determine  the  factor 
of  safety  against  possible  fatigue  failure.  Four-node  linear  tetrahedron  solid 
elements  were  used  to  model  the  hinge  bases  as  shown  in  Figure  43.  ABAQUS 
general  purpose  FE  code  was  used  to  analyze  the  structure  for  the  dynamic  loads 
(reaction  forces)  applied  along  the  two  principal  directions  of  the  base.  Hinge 
forces  along  the  principal  directions  were  obtained  from  the  experimental  results 
for  the  critical  configuration.  All  base  nodes  were  fully  restrained  to  model  the 
fixity  of  the  hinges  to  the  base  plate.  The  maximum  and  minimum  Von-Mises 
stress  on  the  hinges  were  8,000  psi  (55.16  MPa)  and  5,000  psi  (34.475  MPa) 
respectively.  Applying  the  Goodman  line  method,  a  safety  factor  of  7.39  was 
obtained  for  the  critical  gate  frame  hinges. 

Above  computations  also  agree  with  the  AISC  requirements  for  the 
fatigue  allowable  stress  for  members  subjected  to  reversed  loading.  Based  on 
the  AISC  (1989)  strength-design  consideration,  the  maximum  allowable  stress 
range  for  the  wicket  members  subject  to  fatigue  loading  must  not  exceed  22  ksi 
(151.69  MPa).  This  allowable  stress  is  higher  than  the  applied  mean  stresses 
found  for  both  cases  presented  above.  The  allowable  fatigue  stress  range  is 
estimated  for  the  loading  condition  1  (number  of  loading  cycles  in  the  ranges 
20,000  to  1 00,000),  and  the  type  and  location  of  member  connections  as  in  the 
hinge  base. 
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Critical  Response  Simulation 


FE  analysis  of  the  wicket  due  to  the  most  critical  loading  configuration 
was  performed  to  determine  the  adequacy  of  the  proposed  design.  Design 
parameters  obtained  from  a  three-gate  gap  (3GG-L)  flow  configuration  were 
used  in  the  FE  model  to  determine  the  most  critical  response  for  the  anticipated 
operating  conditions  of  the  Olmsted  locks  and  dam.  A  3GG  provided  the  most 
critical  flow  conditions  among  the  proposed  operating  sequences.  Highest 
fluctuation  of  the  wicket  response  and  the  absolute  maximum  of  the  responses 
occurred  at  or  near  the  breaking  of  the  air-pocket  underneath  the  gate  surface. 
For  the  highest  possible  head  difference,  this  critical  condition  existed  at  about 
24  deg  up  from  the  lowered  position  of  the  wicket. 

An  idealized  pressure  distribution  representative  of  the  most  critical 
loading  configurations  was  used  to  evaluate  the  design  performance  of  the 
wicket  by  determining  the  upper  bounds  of  flow-induced  dynamic  wicket 
response.  Experimentally  measured  first  bouncing  mode  for  the  wet  gate 
provided  the  shaft  stiffness  of  the  cylinder  support  of  the  wicket.  Using  a  single 
degree  of  freedom  model  for  the  bouncing  mode,  the  stiffness  for  the  lifting 
system  was  obtained  using  the  equation  shown  below. 

K  =  Mcof  =  ^|^(21.06Ar2;rf  =95,161  Ibf  /  in.  (5) 

In  this  equation,  the  mass,  M,  was  assumed  as  the  experimentally 
measured  equivalent  thrust  exerted  on  the  cylinder  by  the  flowing  water.  The 
added  mass  obtained  from  the  experimental  results,  as  discussed  in  the  previous 
chapter,  was  distributed  on  the  wicket  FE  model  for  properly  simulating  the  wet 
dynamic  behavior  of  the  operational  wicket. 


Idealized  loading 

Upper  boundaries  of  the  random  responses,  based  on  the  Type  - 1 
Extreme  Value  Design  (EVD)  approach  (Wirsching  1993),  provided  the  input 
pressure  loads  for  the  FE  model.  These  estimated  design  values  at  99.99  percent 
confidence  level  were  the  extreme  estimation  of  the  input  loads  for  the  proposed 
flow-boundaries  and  the  likely  service  conditions  anticipated  during  the  life 
cycle.  Experimental  results  for  the  most  critical  flow  configuration  were  used  to 
compute  the  EVD  “design”  values.  Note  that  the  most  critical  response  used  in 
this  evaluation  occurs  momentarily  during  the  raising  and  lowering  of  the  gate  at 
24  deg  position.  The  wicket  has  to  pass  through  this  transition  phase  to  engage 
in  a  prop  support  for  regulating  the  navigational  pool.  A  prop-supported  gate  is 
the  actual  service  condition  for  an  actively  operated  gate.  Service  loads  for  the 
prop-supported  gate  do  not  impose  any  threat  to  the  current  design  as 
demonstrated  in  the  analysis  presented  in  the  previous  chapter. 
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A  linear  spatial  distribution  of  pressure  envelopes  was  considered  to 
represent  the  actual  distribution  of  the  flow-induced  upstream  pressure. 
Experimental  results  indicated  that  such  a  linear  distribution  would  adequately 
represent  actual  loading  situation  for  the  operating  wicket  (Chowdhury,  Hall, 
and  Davis  in  preparation).  The  design  amplitudes  of  the  pressure  envelopes, 
P(t),  were  constructed  by  oscillating  the  random  pressure  response  about  the 
maximum  average,  Pq,  experimental  readings,  such  that  P(t)  =  P^  {t)  . 
Average  maximum  pressure  readings  are  the  maximum  of  the  gauge  readings 
obtained  from  the  critical  flow  configuration.  Figure  49  shows  the  assumed 
pressure  distribution  at  critical  configuration.  A  conservative  estimation  of  the 
downstream  suction  was  obtained  by  uniformly  distributing  the  design  negative 
pressure  (absolute  maximum  gauge  readings)  at  the  downstream  face  of  the 
wicket. 


Figure  49.  Average  maximum  pressure  distribution  at  critical 
configuration  (24  deg) 


The  time  dependency  of  the  pressure  field  was  based  on  a  typical 
normalized  accelerometer-time  history  obtained  from  the  experimental  results. 
An  accelerometer  response  history  included  the  most  sensitive  structural  modes 
and  better  incorporated  the  flow-structural  interaction  due  to  the  flow-induced 
vibration.  For  design  evaluation,  thus,  a  reconstructed  time-dependent  pressure 
history,  Pn(t),  based  on  the  measured  normalized  accelerometer-time  history, 
would  represent  one  of  the  most  dynamically  critical  loading  configuration.  A 
typical  normalized  pressure-time  history  indicative  of  the  dynamic  fluctuation  of 
the  anticipated  pressure  response  was  obtained  from  the  random  accelerometer 
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history,  a{t) ,  using  the  equation,  .  In  which,  p  is  the  mean  of 

V  J 

the  accelerometer  record.  Thus,  the  peak  of  the  normally  distributed  pressure 
signal  would  have  magnitude  exceeding  the  design  magnitude  obtained  from  the 
EVD  estimation.  The  synthesized  signal  would  have  the  frequency  content 
similar  to  that  of  the  original  acceleration  record.  As  a  result  of  such 
idealization,  a  built-in  load  factor  was  implicitly  considered  in  the  critical 
response  evaluation  of  the  wicket. 

Figure  50  presents  a  normalized  accelerometer  time-record  for  the 
physical  wicket  model.  This  record  was  obtained  from  an  accelerometer 
mounted  at  the  top  right  comer  of  the  wicket  during  the  most  critical  flow 
configuration.  For  computation  efficiency,  however,  it  was  required  to  analyze 
the  gate  for  a  fraction  of  the  total  record.  A  2.0-sec  time  record  was  analyzed  in 
the  FE  model  to  simulate  the  dynamic  behavior  of  the  operational  wicket.  As 
shown  in  Figure  50,  this  segment  of  the  record  had  a  maximum  amplitude  of 
1 .33  (ratio).  A  comparison  of  the  auto-spectra  for  the  full  record  (30  sec)  and 
partial  record  of  2.0  sec  indicated  that  the  reduced  time  history  adequately 
preserved  the  signature  pattern  of  the  total  record.  Figure  51  compares  the  auto¬ 
spectra  for  both  time  records.  This  figure  shows  that  the  major  resonant  peaks  in 
the  spectra  for  both  time  records  were  almost  identical  to  each  other. 


Figure  50.  Normalized  random  acceleration  (ratio)  record  for  the  critical 
configuration 
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Figure  5 1 .  Auto  spectra  comparison  for  the  full  and  partial  time  record 


Frequency  analyses 

Frequency  analyses  of  the  gate  at  24  deg  position  modeled  two  different 
spring  stiffnesses  for  the  lifting  system:  95,161  Ibf/in  and  1 1,382  Ibf/in.  In  both 
analyses,  water  was  modeled  by  the  addition  of  masses  distributed  across  the  top 
surface  of  the  gate.  Analytical  mode  shapes  were  similar  for  both  conditions,  and 
are  shown  in  Figure  52.  Calculated  frequencies  are  presented  in  Table  17. 

As  would  be  expected  from  the  results  of  previous  analyses,  the  first 
bending  mode,  mode  2,  and  the  first  and  second  torsional  modes,  modes  3  and  4, 
were  not  affected  by  the  stiffness  of  the  lifting  system.  In  previous  analyses, 
the  frequency  of  mode  2  was  not  affected  by  support  conditions  or  added  mass 
and  the  frequency  of  mode  3  depended  primarily  on  added  mass,  while  the 
bouncing  and  second  bending  modes  were  heavily  dependent  on  support 
conditions.  In  the  analyses  in  the  24-deg  position,  calculated  frequencies  for  the 
bouncing  mode,  mode  1,  and  the  second  bending  mode,  mode  5,  were  lower  for 
the  lower  spring  stiffness. 
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Time-history  analysis 


Since  an  implicit  time-history  analysis  is  computer  resource  intensive, 
the  analysis  was  only  continued  through  2  sec  of  the  synthesized  pressure  history 
record.  This  period  extended  through  the  first  peak  in  the  synthesized 
acceleration  record  and  provided  reasonable  correspondence  with  the  test  modal 
data.  Calculated  accelerations  of  nodes  near  the  hinges,  near  the  center  and  at 
the  top  of  the  gate  are  shown  in  Figure  53.  See  Figure  19  for  node  locations. 
Since  the  initial  0.25  sec  of  each  record  exhibited  large  accelerations  due  to  the 
initial  shock  of  the  dynamic  loading,  this  portion  of  each  record  has  been 
removed.  Displacements  at  the  opposing  top  comers  of  the  gate,  shown  in 
Figure  54  are  the  same  throughout  the  analysis,  indicating  that  torsion  was 
negligible,  and  the  bouncing  mode  dominated  the  response.  Typical  peak 
displacements  occurred  at  2.748  sec  and  2.787  sec,  as  shown  in  Figure  54. 

Figure  55  is  a  displaced  plot  of  the  gate  at  2.787  sec,  with  a  displacement 
magnification  factor  of  20. 

Maximum  principal  stresses  at  the  top  and  bottom  faces  of  the  gate  were 
less  than  12,000  psi,  and  are  shown  in  Figure  56.  Larger  stresses  occurred  at  the 
end  plates  where  the  gate  hinges  were  attached  and  at  vertical  stiffeners  near  the 
center  support.  Maximum  principal  stresses  and  Von-Mises  stresses  at  the  end 
plates  are  shown  in  Figures  57  and  61.  Maximum  principal  stresses  in  vertical 
stiffeners  are  shown  in  Figure  58.  Maximum  principal  stresses  in  the  gate  and 
support  rod  (ball-cup  connection)  hinges  were  also  relatively  high,  with  stresses 
in  the  gate  hinges  approaching  40  ksi  at  isolated  locations  corresponding  the 
regions  of  high  stress  in  the  end  plates.  Maximum  principal  stresses  at  these 
locations  are  shown  in  Figures  59  and  60.  Von-Mises  stresses  at  the  gate  hinges 
are  shown  in  Figure  62.  In  general,  however,  except  at  vary  few  isolated 
localized  spots,  the  upper  boundaries  of  stresses  obtained  from  the  FE  analysis 
do  not  exceed  the  allowable  stress  limit  for  Grade  50  steel.  Since  the  applied 
maximum  principal  stress  of  40  ksi  and  also  the  calculated  Von-Mises  stress  of 
40  ksi  does  not  exceed  the  material  yield  stress  of  50  ksi,  respectively,  the  stress 
state  will  not  cause  any  material  failure  according  to  the  Tresca’s  maximum 
shear  stress  theory  and  the  Von-Mises  maximum  distortion  energy  failure  theory 
(Buchanan  1988). 

Even  though  alarming  magnitudes  of  maximum  principal  stresses 
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occurred  at  isolated  locations,  no  definitive  conclusion  about  their  validity  can 
be  ascertained  at  this  point  without  performing  a  detailed  FE  analysis  of  these 
hot  spots.  Potential  problems  of  element  distortion  and  poor  aspect  ratio  can 
induce  numerical  errors  in  FE  analysis  (Ridlon  1987).  Like  the  element  shape, 
the  relative  element  size  in  relation  to  the  adjacent  elements  can  affect  the 
numerical  conditioning  of  the  stiffness  matrix  (Shephard  1987).  Numerical  ill- 
conditioning  arises  if  the  stiffness  terms  in  a  matrix  differ  by  several  orders  of 
magnitudes.  This  appears  to  be  the  case  for  the  elements  showing  maximum 
stresses  in  the  gate  and  prop  hinges.  As  seen  in  these  figures  (see  Figure  59  and 
62),  the  most  localized  maximum  stresses  occur  at  elements  with  distorted  shape 
and  poor  size  in  relation  to  the  adjacent  elements.  (The  element  aspect  ratio  is  a 
measure  of  the  squareness  of  an  element  and  is  often  measured  by  dividing  the 
longest  element  side  by  the  shortest,  or  in  terms  of  the  angles  of  the  element 
comer).  A  detailed  FE  analysis  by  subdividing  the  base  of  the  hinges  with 
properly  sized  and  shaped  elements  is  recommended  for  determining  the  realism 
of  the  existence  of  high  localized  stress. 
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Figure  52.  Modes  1  through  5,  gate  at  24  deg 


Table  17 

Frequencies  for  Modes  1  through  5 

Mode 

Frequency  (Hz) 

Spring  Stiffness  95,161  Ibf/in 

Spring  Stifftiess  11,382  Ibf/in 

1  (rigid  body) 

37.67 

25.48 

2  (bending) 

57.74 

57.74 

3  (torsion) 

73.75 

73.51 

4  (torsion) 

84.00 

83.75 

5  (bending) 

167.05 

157.91 
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Figure  53.  Accelerations  at  nodes  near  hinges,  at  center  of  gate,  and  near  the  top 


of  the  gate 


Figure  54.  Displacements  at  top  comer  nodes 


Figure  55.  Displaced  shape  at  2.787  sec,  magnification  factor  20 
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Figure  56.  Maximum  principal  stress,  gate  top  and  bottom  surfaces,  t  =  2.787  sec 
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Figure  61.  Von-Mises  Stresses  at  end  plate 


Figure  62.  Von-Mises  stresses  at  gate  hinges 
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Results  and  Conclusions 

An  analytical  FE  model  was  successfully  used  to  predict  the  static  and 
the  flow-induced  vibrational  response  of  the  Olmsted  wicket  model.  This  study 
effectively  identified  the  most  sensitive  parameters  and  the  upper  boundaries  of 
the  responses  for  the  anticipated  operating  conditions.  Boundary  conditions 
were  the  most  sensitive  parameters  among  the  variables  studied  in  this 
investigation.  Added-mass  greatly  influenced  the  dynamic  characteristic  of  the 
wet  wicket.  This  study  also  found: 

a.  The  first  bending  and  first  torsional  modes  were  independent  of  the 
boundary  conditions  at  the  supports,  and  the  first  bending  mode  was  also 
independent  of  head  of  water.  Therefore,  these  modes  were  well  modeled 
in  the  response  spectrum  analyses.  However,  the  rigid  body  mode  and  the 
second  bending  and  torsional  modes  were  heavily  dependent  on  boundary 
conditions  and  added  mass  and  were  less  well  represented  by  the  two 
support  conditions  modeled  at  the  prop  rod  end  in  the  analyses:  full  z 
restraint  and  no  z  restraint.  The  use  of  a  weighted  mass  also  contributed 
to  differences  from  test  data  for  these  modes.  More  accurate  estimations 
of  eigenvalues  can  be  obtained  by  using  the  correct  added  mass  for  each 
mode. 

b.  A  MAC  relation  between  the  corresponding  analytical  and  experimental 
mode  shapes  confirmed  the  similarity  of  the  first  four  fundamental  modes 
of  the  dry  wicket.  A  high  degree  of  correlation  existed  between  the 
modes  compared  for  the  dry  wicket.  The  participation  factor  obtained 
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from  the  analytical  results  indicated  that  the  most  influential  mode 
( mode  #  5)  that  contributed  significantly  to  the  dynamics  of  the 
dry  gate  had  a  substantially  high  correlation  (a  MAC  of  0.984) 
between  the  experimental  and  analytical  shapes. 

A  good  modal  correlation  of  the  corresponding  most  contributory 
modes  indicated  a  good  analytical  representation  of  the  wet 
physical  model.  The  strong  modes  that  are  most  contributory  to 
the  dynamics  of  the  wet  wicket  had  a  high  degree  of  correlation 
(a  MAC  of  0.826)  with  that  of  the  physical  mode. 

c.  The  static  test  was  modeled  fairly  accurately  using  a  spring  stiffness 
of  15,1761.6  Ibf/in,  although  the  predicted  displacements  at  the  top  of 
the  gate  were  22.8  percent  less  than  test  displacements,  indicating  that  the 
F.E.  model  may  be  slightly  stiffer  than  the  test  1 :5  model. 

d  For  the  time  period  analyzed  in  the  65-deg  position  of  5  sec,  the 

frequency  content  of  the  dynamic  analysis  acceleration  records  differs 
somewhat  from  that  of  the  test  accelerometer  records.  Amplitudes  of 
structural  frequencies  were  low  due  to  the  low  level  of  excitation,  but 
peaks  are  visible  in  all  records  at  approximately  46  Fiz,  55  Hz,  and  72  to 
75  Hz  that  may  correspond  to  the  rigid  body  mode,  the  first  bending 
mode,  and  the  first  torsional  mode,  respectively.  ABAQUS  amplitudes 
are  generally  lower  than  test  amplitudes  for  the  first  and  second  mode  and 
higher  than  test  amplitudes  for  the  first  torsional  mode. 

e.  Calculated  reactions  (and  therefore  stresses)  are  extremely  dependent 
on  boundary  conditions  at  the  prop-rod  end,  and  test  reactions  could 
not  be  reproduced  by  a  single  ABAQUS  analysis.  However,  the  two 
time-domain  dynamic  analyses  (z-restraint  and  partial  z-restraint  at  the 
prop  end)  produced  reactions  that  bounded  test  reactions.  In  both 
analyses  for  the  prop-supported  wicket  maximum  principal  stresses  were 
low  and  should  not  present  a  problem  for  the  simulated  operating 
condition. 
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A  critical  dynamic  response  analysis  of  the  FE  model  indicated  that  a 
maximum  principal  stress  of  about  40  ksi  would  develop  on  a  very  few  isolated 
localized  spots  on  the  wicket  for  the  extreme  loading  condition  expected  for  the 
proposed  flow-boundaries  and  the  likely  service  conditions  of  the  Olmsted 
wicket.  In  general,  however,  except  at  these  few  isolated  spots,  the  upper 
boundaries  of  stresses  of  about  12  ksi  on  the  top  and  bottom  faces  of  the 
wicket  do  not  exceed  the  allowable  stress  limit  for  Grade  50  steel.  Based 
on  the  Tresca’s  maximum  shear  stress  and  Von  Mises  maximum  distortion 
energy  failure  theories,  no  failure  state  is  anticipated  due  to  the  simulated 
applied  load. 

Larger  stresses  occurred  at  the  end  plates  where  the  gate  and  prop- 
rod  support  hinges  were  attached  and  at  vertical  stiffeners  near  the  center 
support.  These  hot  spots  are  believed  to  be  due  to  numerical  error 
resulting  from  a  few  elements  with  distorted  shape,  poor  aspect  ratio,  and 
relatively  smaller  size  element  compared  to  the  adjacent  ones.  Like  the 
element  shape,  the  relative  element  size  in  relation  to  the  adjacent 
elements  can  affect  the  numerical  conditioning  of  the  stiffness  matrix 
(Shephard  1987).  Numerical  ill-conditioning  could  be  a  problem  if  the 
stiffriess  terms  in  a  matrix  differ  by  several  orders  of  magnitudes. 

A  detailed  FE  analysis  by  subdividing  the  base  of  the  hinges  with  properly 
sized  and  shaped  elements  is  recommended  for  determining  the  realism  of 
the  existence  of  high  localized  stress. 

This  critical  analysis  used  the  upper  boundaries  of  the  pressure 
responses  based  on  the  Type  - 1  EVD  approach.  These  design  loads  at 
99.99  percent  confidence  level  provided  the  most  critical  estimation  of  the 
input  pressures  which  would  occur  only  during  the  lowering  and  raising  of 
the  wicket  at  or  near  the  separation  of  the  air-gap  beneath  the  wicket.  This 
analysis  also  included  a  pseudo  load-factor  of  1 . 1 8  in  the  analyzed 
segment  of  the  input-pressure  time-history.  The  wicket  has  to  pass 
through  this  transition  phase  to  engage  in  a  prop  support  for  regulating  the 
navigational  pool.  A  prop-supported  gate  is  the  actual  service  condition 
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for  an  actively  operated  gate.  Service  loads  for  the  prop-supported  gate 
would  induce  a  maximiun  principal  stress  of  2.5  ksi. 

In  general,  the  predicted  critical  stress  for  the  Olmsted  wicket  is  far 
below  the  AISC  design  flexural  stress  of  33  ksi  for  the  compact  I-beam. 
This  comparison  effectively  validates  the  design  of  the  Olmsted  wicket.  A 
fatigue  analysis  using  the  modified  Goodman  line  method  also  indicated  a 
very  high  safety  factor  of  7.39  for  the  critical  gate  frame  hinges.  This 
observation  is,  however,  limited  to  the  operational  sequences,  head 
differences,  and  geometric  configurations  of  the  wicket  and  the  flow-field 
covered  in  the  experimental  study. 
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